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SYl'JOPSIS 


Ki®.n; G-m® iossfa 
M.Tech.(Mecb.) 

Indian Institute ox Technology Ilanpur 
3eptei±ier, 1973 

INSTITUTE LMTSL PLillMIIK FOE .ISMISSIOHS- 
A QU/ATTITATIVE APPROACH 

The 77ork presented here pertains to an educational 
planning proUlen of admissions for postgraduate students for 
Master of Tech.nolog 3 !’ at Indian Institute of Technology, Kanpur. 

Tv/o models hare been deve lope d^or foreca,sting the nmtber of students 
who E,re likely to join the various engineering disciplines. The 
proposed models are based on the Linear Cyclic Forecaster and 
Ealman Filter Algoritlm, In the linear Cyclic Forecaster, the 
nuidoer of students who are likely to join the various disciplines 
is considered as deterministic while in SoAnan-Filter it is 
considered as probabilistic, 

Dynaric progrrmnii^ has been xised for determining the 
number of seats to be provided for the various di3cipl5.nas during 
the forthcoming periods or semesters. The '■■■ntimization 7;as carried 
out using the criterioix of mnimizing the cost function subject to 
various constraints. The constraints considered are : Financial 
resources available and the pemissible. values of the mnimun and 
maxii-jum numljer of seats available' in each discipline, 

u The proposed is have 'been validated by usiiig^ the 
data; furnished'bydtite;::Ppst~grhdaate. Office ; 0 Eanpur* ;/ 


, ty-T/ tD- 


IlITBCSUGT'IOn 


1,1 CSIIISSAL IIITEODUCTIOIT 

It is mdely accepted that some neasrire of echicational 
planning is necessa.r^^ to ensrxe that the edncatioiaal system discharges 
its hea-vy responsibilities to the society in an efficient mamier. In 
India, the expenditure on education ferns a significant portion 
(a-boiit S^) of nation*s economic re source s» Further it is believed 
tha,t if planned properly, education mil sake very important contri- 
butions to the econoEiic gro^fth and the social development^ of the 
country# Due to the complex nature of our society, there are very many 
conflicting objectives which the planner has to bear in xaind# These 
conflicting objectives nay result from i) rapidly increasing demand by 
the society for literacy, ii) equality of opportunity, iii) laanpOTOr 
reqydremerits for the growth of economy an-d iv) individual's freedom 
of choice in the type of education he/she wants to pursue, . The belief 
that education is an important contributor to the , economic growth 
ha.s very recently triggered -a lot of , interest amongst researchers to ■ 
come, up with methodologies for , effective’ educational' planning* 

: ■ Soj^ of , the ' important 'Vquestions ' relat^^ to the, educational ,' 

planning r'are;' h, f.' 

i") ' lifhat' amount 0ffs.oeiety'*s ' resor^^^ /should' , be 


to education? 
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ii) ?Jhat should be the ratio amongst the educated people 

of various disciplines keeping in view the requirements 
of the country? 

iii) Eow should the total resources be distributed among 
various types and disciplines of education? 

iv) Ifhat would be the likely rate of return for the expen- 
dittire incurred on an educational system? 

These questions pertain to the economic efficiency of 
educational planning of the country as a miole. Since educational 
systems are very complex in nature, no attempt has been made in the 
^aast to find optimal solutions for education planning problems on 
national level. Keeping in view the intricacies of educational 
planning, many research workers have tried to solve some of these 
problems at institute or university level by various techniques 
like simulation \20j , ?EST/CM f^28j, linear programming f 21, control 
theory [23, 3^ , etc. These techniques have primarily been used 
to study the flow of students through an educational system without 
considering any resource constraints, allocation of resources to the 
various facilities of an educational institution, scheduling students 
to classes, development of curriculum, etc. 

However, it seems that no— body so far has studied the 
apportioning of financial resources keeping in view the interaction 
of various disciplines among which these resources have to be 
distributed. Further no research has been reported for forecasting 
the requirements of the various disciplines taking into account the 
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effect of : 

\ ^ 

i) salary gradient for the graduate of various 
disciplines 

ii) the number of students called for interview in 
these disciplines 

In this thesis an attempt has been made to study a sub- 
system of the total educational planning system. Ihe subsystem 
selected is the planning of the post-graduate technical education 
at the Indian Institute of Technology, Kanpur taking into account 
the effect of above mentioned two factors. Mathematical models 
have been proposed and tested to tackle the following problems : 

1. Forecasting for enrolment in the various disciplines 

• of engineering* 

2, Number of seats to he provided in the various disciplines 
t6iking into account various constraints. 

For forecasting two models have been proposed. The first 
one is a linear cyclic forecaster and other is based on the 
principle of Kalman Filter. For predicting the number of seats in the 
various disciplines, dynamic progranning is used to advantage^ 


This is the difference between the average salaries offered to 
a po stgraduate and an undergraduate student. 



CHAPTER II 


LITEBATUBE REVIEW 

The review of the literature on the educational planning 
suggests that the models can he grouped into the following four 
categories, 

2.1 CONTROL THEORY TYiE MGBEL 

Koenig and Keeney f_23} developed a model which can he used 
at the university level for the distribution of educational 
resources to efficiently seive a number of students, Alper and 
Smith [33 developed a model for finding number of places to be 
made available in a sector of education to satisfy an unknown 
social demand. Besides the control theory concepts, dynamic 
programming was used to determine the seq[uence of decisions which 
optimizes the cost criterion. The major assumption in their work 
is that the number of potential entrants is limited to an unltnown 
constant level. One might question the feasibility of such an 
assumption, 

2.2 SCHEDULING MODEL 

Blackesley {123 developed models for scheduling students 
to classes at the beginning of each term, in data processing of 
student grades, records, etc. Abell Ll] developed optimum seeking 
algorithms to schedule final examinations of an educational 
institute, Taft and Eeisem L37j developed a mode 1 which deter- 
mines the sequences of courses, for a student in an academic^^^^ 
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institution which naxiudzes mstery levels at graduation* 

2*3 ig)DEL 

Myriek, developed a network model of the Piirdue 

Freshsssn Engineering curriculum* The model has heen programmed 
for simlation using the GASP II procedure* Curricultua, , f or a 
student can he considered as a time based precedence graph, v;here 
the arcs represent the event the-t a course has been completed. 
Figure 2-1, for example, represents one of the t'tvo semester 
curricula available to a student in the Purdue Department of 
Preshoan. Engineeringv A‘ network.. .cac--repre'sent ''ah'" individual 
'stiident or a group of students-~talEing the same curriculum* 
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2*4 OTESil imims 

Ju&y 1201 developed simalation model to solve educational 
problems of resource allocation* Sanderson T3G1 developed a model 
to study the impact of increased enrolments on 'educational faci-* 
lities* In 1967, educators, statist! ciani}, engineers c-nd other 
professionals met in Washington in the first 'IJ*S*A* conference 
devoted e^iclusively to the application of operations research 
tecimiques to the field of education* 

Pritsker £303 used decision theorjr approach for enrolisent 
policies of iirisona State University* His objective was to predict 
the number of students expected to be enrolled in future* For 
this he minimized the total cost keeping in view a fixed financial 
budget* Different penalties were specified for poor predictions - 
both overprediction and underprediction* The two terms used abcve 
and the associated penalties in- terms of cost are delineated in 
the subsequent sections* 

1 • Undernre di ction : ■ 

The term underprediction implies a predicted value less 
than the actual value* It ?/as further divided into the , following 
twO' categories ' 

a) ' Major , uiiderpredicti on :■ If "the, predicted, value is "less than 
a specified percent of actual value, Pritsker defines it ,as major 
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underprediction, 3ero4n a request has to be mde to grant funds 
ovei" and above the schedule! anount# Thus a cost/iinderpie dieted 
student uas indicated by a constant C. , 

b) I/Iinor unde rpre .diet ion : If the predicted value is less than 
a, specified percent of actuad value but greciter than raajor under- 
prediction, Pritsker defines it as minor undGrprediction. Herein 
no request has to be loade to grant funds over and above the scheduled 
ainouat and a penalty?- due to opera-tion on a lowr budget is assessed# 
The cost/underpredicted student was indicated by a constaat Cg. 

2# OverT3ie diction 

The term ove rpre diction implies a predicted value greater 
ii^ n c mOi o ne actiaal valuer It Tmsflirtber divided into the following 
two ca.tegories : 

a) i/Iajor overprediction : If the predicted value is greater than 
a specified .percent of fictual value, Pritsker defines it as major 
overprediction» Here-ia a request has to he made to refund funds 
over and above the scheduled amount. Thus a cost/overpredicted 
student was indicated by a constant Cg, 

b) I/Iinor overprediction : If the predicted value is greater than 
a specified percent of actual value, but less than major over- 
prediction, Pritsker defines it as minor overpredictipn. Herein 
no request has to' be made to refund funds over and above the 
scheduled araount. Thus a cost/overpredicted student" was 
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indicated by a constant 

An analj^tical a.pproa.cli gave the best predictor of 
enrolment to be made in future# L sensivity 'anaijrsis was 
carried ont for various values of = 0,2, G#4, 0#6., Q#8, 

1#0, 1#5, 2.0, 3.0, 4.0, 5.0, 9.0. 

A recently developed technique of Kalman filter based 
on the principle of control theory can also be used for predic- 
tion of enrolment^ Kalioan (21} considered the linear estiFiator 
for GstiDiating the states of a nonlinear control system. EaliBa^n 
defines noise corrupted observations in a system as a nonlinear 
system. Salman considered’ the linear estimator as a linear 
2 ?unction of observation data and stated that the best linear 
ostirisator is best of all estimators (linear or nonlinear)^ Esti- 
mation of the state in the. presence of unknown parameters of the 
nonliraear sj^-stem ?/as trea.ted hj I&gill£241 and Ililborn and Lainl0tis|183-i 
Cox 118 3 ^>^d Balkrishnan 16 3 treated, sar^ problem via a maximum 


likelihood criterion^ 


CEAPTEB III 


FaOBISM FDEi'ULfi'.TIOI-T 


3.1 IMTEODUCTIOII 

This study on the planning of an edncational aysten for 
a technical institution involves the prediction of the number 
of seats to be provided for the various engineering disciplines 
keeping in view the socio-econoiaic constraints. For this the 
objective function at any period k, is expressed by tiie 
equation : 


. f (jjj , xJ , cj , cS , pj , 3j , G^) 


(3.1) 


Wile re , 




Z =: the value of the objective function in terms of 

iC 

th 

cost units for the k period, 

. , , ■ who are 

« the number of students£likely to join the discipline 

j for the period, 

U? s= the number of seats to' be provided in the discipline 

„k ■ ,V.;, ■ 

th ■■■ 

; j '.f ,or'"the' k period,' ■ 


■^oln', of' 'p^ be ' a' seaie stern' 
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’"2 


= expendittire of the disciplir»e j for the period, 


•R 


n 


s= the penalty in terms of cost units for Timder- 

■X- 

pr edict! on for the discipline j 
= the penalty in terms of cost nnits for over- 


¥ 


prediction for the discipline J 


th 


= the facnlty required for the discipline j for the 


1 -1 
k T3eriod, 


th 


the hiiii dings of ty 3 e required for the k ""period 


= the buildings for the postgraduate students, 


the tota.l nur-iDGr of disciplines ta!^:en, 


IT = the nurber of periods in the planning horizon< 


The nuiiiaer of students called for interview is a function of 
nany variables ♦ These variables are listed below : 

1, The nurfber of applications received 
2 ^ Academic perferaBnee of the : students in the previous 
enmei nations 

S# ilcadeiiic standards of the various institutions from 
77hich the applicaoffcs have graduated 
4#. The ntuaber of seats :available., for the semester for Tvhich 
' the forecast is: to, b^e. ■ imde..keevping''in view' the 
■„:of' student S', who would’’ 'be .gradiiat ing tte start 'of this 


■*'' 'T'h'i's'' " occurs' when^ "the hurb'er ■ofTvs'tudents' '^who ' are; likely , 'to': join 
, a ,, di,s'cinliiae' '' are' th'aii': tho'' ’.Biinber,:,'.of ■ 'sect.,s vprovidedx, 'f or ' ,■ 

, '';;,'that,'''discipiine#v , ; 

■■^*Thi,Sv;'occurS' when' ''the'' 'uiiiiber^ are,',", like,ly' to, ; .Join ;; 

'a,',disci'nli,ne ''■ bo,,^''' 'than^ ' the;.\,nu^ scats, provided 'for" that 

discipline# 
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5, Total hostel accouimodation available. In case of 
completely residential institutions, the number of 
hostel seats vacant is a very important factor for 
consideration. 

It is not possible to develop an explicit mathematical 
relationship for the exact prediction of the number of students 
•(vho would be called for interview because it is a function of 
many intangible factors (factors 2 and 3 as listed above). The 
difference between the nunber of students to be called for inter- 
view as found hy using the mathematical model and the actual number 
of students called for interview is called error in prediction, 

Por all the practical purposes one can assume that this error in 
prediction follows a nonaal distribution because it is the net 
result of very majtiy independent intangible factors. 

The nuEiber of students who are likely to join in the 
various disciplines for any planning period depends upon the nurijer 
of students called for interview in the past periods, the nvoniber of 
students joined in the past periods, salary gradient for various 
disciplines during the past periods and the reputation of the 
institution. The difference between actual and predicted value 
is the error in the prediction and it is assumed that it follows 
nomal distilbution. 

The objective function expressed by equation (3,l) is to 
he optimized subject to the following constraints I 
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(i) Eeso’j-rces available to the institution 


'^k 

j=l 

T/hsre , 


( 3 . 2 ) 


€, = total funds avail 0 l 3 le to the institution in the 

K 


th 


oerioa# 


(2) The total nuaber of seats* sanctioned to the institution for 


the neriod 


^4 -h 

j=i 


(3.3) 


T7b.ere , 


"k 


total number of seats sanctioned to the institution 
for the period* 


(3) The muiber of seats in each discipline at an^ period k should 
fe>il het^Teen its mninaim and EiaxiEsiim values* 

4 ^ 4 -t -i < 3 -^) 


wnere . 




' and are the mnimm and imximin value on the 


th 


number of seats in. the discipline j for, the,, k ; 'p^ 


lPhe,''''^'t ' :seats :tsanctioBed,' ,t'o 'the^'institu^ any';:' 

period' h, is 'E'' function ■ of'.' the total'*' funds; ;schol,ar,sIilp ,' 

to 'the, students and/hostel, facilltie.S'.*',''^,'' 
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(4) The ndniisii!!! and Baximn nnmher of peri'iiGciible seats in eacla 
discipline for four conseentive T3eriGds is bounded* 


E^. + u,^ , + "d o ■<- -d 

Eaa ' k k+1 k+2 k-; 


TTi J 


(3.5) 


where, 

. and are the ninioun and naximia nunioer 01 

imn m.x 

peririssible seats in discipline j for four consecutive 
periods# 


(5) 'Additional constraints : 

(a) Constraint due to buildings - Buildings of type i 
required in period k (S^) less than or equa.1 to 
initial stock plus new construction since beginning of 
planning period k# 

(b) Constraint due to faculty strength - Teachers with 
speci^ility in discipline j' required to accoimodate 
the enrolments in period k (T’^) for discipline. ’J’ 
are less than or equal to teachers surviving' f rom the 
initial strength plus new teachers recruited during 
the planning period# 

all these, restrictions are simltaneously , taken into , 
account ,. the" .problem' becomes ■ too complex to .handle#' In /order to...' 
simplify the problem 'it is /Bece''SBary. tO'/malie . a' few,, assumptions*., 'Tn 
the present woiA% .. only '.the .'constraints'' ''2, .,'3, 'and: 4' „are ' co'iiside red 
,The entire . problem has 'been split into ''two, phases * 
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The first phase handles the forecasting of the nanber of 
students who are likely to join a discipline in a given period. In 
the second phase the optimn noriber of seats to be ’provided in the 
various disciplines is determined taking into account the various 
constraints. 


For the first phase two types of forecasting models have been 
developed. In the first model called linear cyclic forecaster, 

(the number of students who are likely to join the discipline) is 
taken as deterministic quantity. In the second model, is treated 
as probabilistic in nature. This model is based on the control theory 
concept. The details of this model are given in following section. 

In the chapter of Results and Discussions the results of both the 
models have been compared. For the second phase dynamic programming 
has been used for the determination of the optinnm number of seats to 

be provided in the various disciplines, 

3,2 FOBSCASTING THE NUMBER OF STUDENTS HKELY TO JOIN 
3.2,1 MODEL I : Linear Cyclic Forecaster 

This model has been developed to forecast the number of 
students who are likely to join the various disciplines. In this 
model the mriber of students who are likely to join the various 
disciplines is considered as deterministic . The model is developed 
as a tine aeries and is capable of generating demand trends 
a wide variety of slopes, growth rates, cycles and variations. 
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T;:ie tine series used in the Eiodel for generating the 

forecasted demand for the mdber of students mho are likely to 

* 

join the marions disciplines is e:rpressed hy the ecaiation 


JD(t) = A -!■ 3t -}- D(Cos ( -nt/E))+F(Cos{-n t/s) )+H(Cos(T5t/U) )+3S¥ 


T/here , 


JD(t) 


A 


3 

D, P,H 

E, G-,U 

mr 


(3,6) 

K generated desBaiid of the nunber of students who 
joined the discipline for period ’*t”, 

=: the constant which represents the deraand at time 
period zero (”y*’ axis intercept) 
ss Slope of ^’Straight line”, 

=: Coefficients Tjiiich detemine the anplitudes of 

first, second and third cyclic trends respectively, 
= Coefficients which detemine the periods of first, 
second and. third cyclic trends respectively, 
s= a randomly chosen deviate from a norml popuiation 
having a nnan of zero a-nd a standard deviation of 
:SIGLiA. . This normal ■ distribution is truncated at 
three , standard deviations below ■and above; the mean#^ 
' = The nunber; of; available past demand'data, , points:# ^ 


The',, 'above :;equation^ to each of the engineering 

disciplines and tliergfore the superscript j on the variables 
has not been put# 
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The deaaEsd is generated asstiiaing a linear cyclic growth. 

The nanber of cycles occuring in the y.ast data for each 

discipline is detenained by using the technique of spectra.1 analysis ISl} 
and it was found that for each discipline a rsaximna of three types of 
cycles can be inpcsed on the deajand trend. The arbitrary conponent of 
the tine series {D^TJ) is calculated using iZonte Carlo Simlation 
Technique \,25b . This keeps the usual off control points in the 
frec'UGLicy ansilysis of the deriind pattern at the confidence level of 
9G.737S. Using this model a, value for each parameter (A, 3, C, ,,, etc^) 
can be computed for each discipline. In the present work, it is assumed 
that sufficient past data for the number of students who joined in the 
various disciplines is available, and therefore regression analysis 
is used for estimation of parameters. This assumption is validated in 
Chapter IV. 

3 .2.1.1 Solution procedure 

The va.rious steps involved are as follows ; 

STBI X !' An iterative scheo's has been employed to evaluate the 
values S, G and U, In this scheme one takes various 
values of S, G and U in the range from 1 to H/2 in •' 
steps of 1/2. . This is done because the periodicity of 
c;?'C;les have integer values and the minimum possible 
number of periods in: a cycle is- 2, Total absolute error 
in the f oreeaBting in, the- past N- periods is used a,s 
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3TE? II : Once the periods of cycles are identified, the values 
of other coefficients are calculated in a nanner given 
beioo; : 

i) lie defining the variables of tire series to reduce the, 
co:::niita,t iona.1 efforts, as follovs : 


^There 


or 


lr(t) = - DMd 

k'(t) = - mm . (3.7) 

h = t - n/2 


ID{t.) = actual denand for period t 
DIS'G » average of past demands over last N periods 

ii) After substituting equation (3.7) in equation (3.6;, 

the values of the coefficients in equation (3.6) are 


deternined by solving the set of sioiltoiiieous equations. 
The set of .ginultaneous equations is solved by Pivotal 
Condensation Uethod 728] and the values of the elemental 


dHewnlnonU 
la at- r - io es Bl, 


D1 , 31 , 


FI , G1 and HI are calculated. 


The solution procedure is given in Appendis A. The 
technique of Pivotal Conciensa-tion li'.etnod fina.lly 
results in the equation : 


k»(t)*31 - D1 + h * 31 - S'! ^ Cos 


-rrt 

E 


+ Gl* Cos 


TTt 


-a*cos 


TTt 

U 


k • (t )= (Dl-Sl*h+Fi# Cos ^ - Gl*Gos 


G 


(3 #8') 

Hl*Gos . 
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Since k’(t) = JD(t) - Dl-TS- (iby definition) 


or 


or 


h = t ~ n/2 (by definition) 

<]l;‘(t)-DiiFG = (Bl-El*h+Pl Cos 

JL{t) = If- + 


G1 Cos ^ +I-I1 Cos ^)/B1 

*jr U 


SI /, I-K Pi „ TTt 
B1 ~ 2''^ B1 E 


G1 _ TTt 
21 G 


El 


B1 U 

On cenyaring tloe coefficients of equations (3.6) and (3,S), one gets the 
vaines of tbo donsnd ■ocraneters as follows : 


TTt 


(3.9) 


A = (Di + El * N/2)/Bi + Dim 
B = -El/Bl 
5 . M/31 

1? = -Gl/Bl 

H = Hl/El 

STE? Ill : Elie value of the standard deviation of deraand, is 


cs,lculated using the equation^ 


Ali, standard deviation =; 


* (Error i n f o re c as t ) ^ 


IT-f 


, when IT > f 


where , 

f = degrees of f re edon (5 in this co.se) 

and Error in forecast = difference betxTeon actual and 


f ore caste d demnd. . 

Since the degree 'of freedon is 5, tberefore, to use the 
above runtioned nethod a laininiaEi of 6 periods^ denand history raist ;be 
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knovm. As the nnnber of past denand data increases the accuracy 
in the estination of denand paraneters increases, 

STEP W : The values of E, G and U 'which are calculated hy . the 

iterative scheiae are checked by correlation and spectral 
analysis fsij , These techniques are being used to 
calculate the approximate value of periods of cyclic trends 
for situations when one is faced with a large anount of past 
data. The iterative schene as oentioned above becones 
too curibersone to use under such circunstances. 

STEP V : During each period the forecasted demand is compared with 
the actual demand to evaluate the effectiveness of the 
forecaster. Higher the value of the error j the lower is 
the effectiveness of the forecaster. 

Having computed all the coefficients of the time series, 
the forecaster can be used for planning parposes, 

3.2.2 MODEL II : Kalman Filter : 

This model is based on control theory concepts of Kalman 
Filter and t alee s into account the probabilistic nature of the number 
of students who will be called for interview in the various disciplines. 

Taking mmber of students called for interview as the obser- 
vation and the nuraber of students who are like ly to join the 
discipline as the state of the system (mdel), the Kalaan-Filter ^ 
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model can be expressed matlierjaticaily as 


(3.11) 


* "Ic 


(3.12) 


•r "DIa / \ 

= the observation at the k "" period ^Input to the sodel), 

= the state of the syoteri at the k period (probabilistic 
qnajatity)* It io an m coaponent vector (output of 
the model), 

th 

= the detemiinistic salo.ry gradieiat for the k , period 
to the model), 

th. 

1 /,^ = the noise in the state cf the system at the k period, 

th 

» the noise in the observation at the k ' period, 

th 

g. = the coefficient of the salary gradient tern at the k' 
■period, 

h = the coefficient re lasting observation and the last 

th ' 

component of the state vector at the k ^ period, 

th th 

■ dp,, V- '■ ^ . component of the state vector at the k " period, 

in,];!' 

:F . 5 = an mm coefficient UTatris: relating t?ro consecutive, - 

states# , 

Sack of the .terims'''' , 3^ 2 ^,- and . in em (3»ll) and, (3*12) ' ■■■ 

: ' ^'-th ' 

:is a m-dinensional vector* The current, system' state for 'k-fl ; ■ period'" ' 

” ' ' th. ' 

will be mostly affected by the system, state; for k period sad as tbe , 


Tiie E'3ove set of eaftiatioEis { 3 »ll) and (3.12) are applicable to each 
of the disciplines and therefore the superscript’ j on the variables 
has not been put. 
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age of the date for 3 past periods increases, its influence on the 
S 3 ?'stem state for cTirrent period decreases# Therefore, it lias ’been 
assiiined. that the systein state of only past periods ‘affects the 
current state of the systen. For the sake of simplicity considering 
only ^-dimensional vector case, the equations (S.ll) and (3*12) can 
he Tjritten as : 


J 

’ 0 

1 i [ X , i 
\ l»k i 


^ “S ! 

* 1 i 

: 1 ,k 

II 

1^1 

\ i 

f 1 i y ' 

2 } L 2,k . 

i ■ 

L 

\ 

1 

L i " 

• ^ J 

' V 

1 2,k j 


(3.13) 




?,,k k 


(3.14) 


Tj-hei-G , 


X. , = tiie first coaDonent of tke state vector at the 

1, k k 

th 

k " period, 

X = the second componeBt of the state vector X, at the 

2, k K 

th . , 
k period, 

f . , f„= eleneutc of the coefficiert Batrix F, 

Souations (3,13) and (3.14) constita.te the fundaiaental equat ions for 
each discipline of the edncatioBal systeia under consideration. 

Given the observations -g, "^2’*'"’^ deteimnation 

of the estimate of the systeE state for say period K^, 0 ■$ K, 

the prohleia ox predicting the S3ratem state is called "dc.ta sEoothing" 
nrohieEi* For Kj = k and k, the prohlons are referred as: ■ 

."filtering" aha "forec.asti.ng" prohieEas respectively. 
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Defining j 


= p (y“, g, 


{3.15) 


pCif") p (X,, X„, ... X 


^ P ^ 1 > ^“2,1’"'*’ ^"1 , n ’ ^"2 ,.n 


T^e "T-xiizdsaticn io "^erfonx-d over it. , ^ i = 1,2, 

i,K 

f » F ii» Cn Gimlifying equation (3.15)', the lollomng 
1 2 

rocursive rGlc^tionsIaip can be obtained : 


\ = P V-f„k) pVX- 1' ‘l’ ^2’ ®- X-1 


(3,16) 


The ref ore 

log fu 

notion of for 

n 

tically : 




In L s: *-»C 

n 

n 

Z 

0=1 

(Y.-h .f 

,1 2,.1 

n 

0=1 

■ ■ 


:, (x„ -f, X, /. -f„ x„ . .-gs . )‘ 


o 2 

2 <r 

9 r 


f Tr -^r \ ^ 

XT o 2 

j=l 2 cr 


(3.17) 


T?!lQrC , 


=s Ta.riaB,ce of the error , in the observation, , 

■ '2.: ":2'": 

cr- CT" - ElecffiEts of tbe v:;riai3CG in the system state 

h’ ^"2 


'G, =: { ' constant ' tein# 
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The objective is to find oat the systen states and 
paraneters. For obtaining naxinon likelihood estimate of the system 
states and values of the paramgters f^, , fg, g and h, the joint proba- 
bility ftinction of Xj^ and will be maximized. For this equation(3.17) 
is differentiated v/ith respect to X. , , i=l,2, f.,fn,g 

and h. Each expression is separately equated to zero. This results in 
a system of non-linear equations in the variables X. , , i=l,2, 
k=l^...,n, fgi g This system of non linear equations can 

be separated into two sets of equations, that are linear in corres** 
ponding unknowns, namely, the states and the pa razoeters. 

Performing the maximization of equation (3,17) one gets, 


(in L )=f, 

cih.k ° 1 


(X 


-f.X. ,-f.X. 


-gS, 


2.k+i -ri,k -r-2,k 
Ov? 


\ - 


<Jv? 


=0 (3.18) 


k = 1 ^ # # jjQ"*! 






■2,k 


(in L ) 
' n' 


« h 


Y -bX„ , 

n 2,k 


^2 
^ W 




^.k~^l^l .k-l~^2^2.k-r^\ -l 




^2 .k+l~^l^l .k'‘^2^2.k'“^\ 


(fv. 


1 fk+l 2 ,k rt 

(Tv 

1 


^^,n 


(in Ii ) 
' n' 


X^ - X„ . 
l.n 2.nr-l ^ 0 

uv^ 


(3,19) 


(3.20) 
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Y-Ii Z. 


. -f. Z. 


-f JZ 


i (in - 2,n 1 -j ,n-l 2 '^2,a -l ^ n-1 ^ q (3^ 

’•c ^ ^ 2 2 
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-1— f,. X 1 V 

t>-i "J" Zj 


j=l 


n 


^l,.i- l~ ^2^'2..1-1~^'^.1-1^ 
2 

rVo 


z. ._^=o (3.22) 


L.(in n 2, » ,1-1 2^. , .1-l~g^ j-1^ 


'3 e 


j=i 


^ h ■^J ■ 


(Y-h .) 
■IJ— -iLtil 

2 


2 

(Tv, 


X, . = e 
2,3 


j-r 


(3.23) 


(3.24) 


j=i (f^ 

Tlia sot of eciur.ticns (S.lo) to (3.24) can be separated into t^ro sets 
of e-qr. ations that -are linear in corresponding ^^n!'no’l'ms , nanely, the 
states and pOiraEctors f^, g and h. 

The set of. equations (S.IS) to (3.2l) can be rewritten in t'se 
inatrix fora as 


3 


n n n 


(3,25) 


whore. 


.Jx J 


El r I 

\ ’ h> J 











?,6 









^ 1 ^ 1,0 
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The set of. Gcpaations (3.22) to (3.24) can be rewritten in natrix 

f orn as 


■5 P = r 
n n 


T7.IQrG' 


P 


(^2V 


Ch 


(3.2S) 
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It can be noticed that the set of ernaticns (3,25) are 

linear in Z. ^ and can be solved to give the 

v-alnes of provided the values of g and li appearing on 

the right h^n:id side of the equation (3*25) are knovn. Sinilarly 

the set of aquations (3#28) are linear in f ^ g and, h and 

requires the vaMes of Z. , (i = 1,2, k = in, order to be 

X u 

solved for f ^ , fg^ g and h* Therefore, these equations can best., be 
solved iteratively using a suitable nurierical technique* 

The airplication of the principle of Kalnan filter to the 
educc..t iona. 1 planning problen has resulted in values of the system 
states for future ’:)©riods# In the a]3sence of a,ny constraint the 
values of the sj^sten state for future .periods could be taken as the, 
nuribei* of seats to be provided at thsi-t period* IIoTjever, , the cons- 
traints imposed on tlLO problem iiay not alio?/ to take tne values of . 
system stcites a.s the number of seats to oe provided. ,, Theref ore , a 
suitd:.;.lG stage by , stage optimization technique is reqxiired to yield 
the final results. ' In- the present circumstaHces the dyn omie ;pro.gra-. 
mruBg seems to be a powerful optimisation technique to ,geo 

optimal , decisions • , ■ ■ . , , , 

The .results ,of /Linear ■■Cyclxc 'Forecaster and Zaimn Filter , 

' te chBique are ■ given . 'in,. Chapter b' . d : ^ 
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o#c^ ,j.,..L':Jr.^-..u ^ .. ..^^- ;:.‘^.ij.i.i.ii.^ i.HU.:J:::ijjLi 

Siiios t’:ie T’rediction of the opt iniin nuiiibe r of scats to 
he ;:::c'/iaof ia the various disciplines for each period is a stage 
by stage- p*rcbie::: c:-: prediction subjected to Taarioois const ro-ints, 
theTOfoT<& dynaiaic prograrieing has been selected as a riathenatical 
tool for to-hing saca decisions* 

in tae* present disserta.tion the ontiinization is cariried out 
nsinp t::o criterion of rniniiiiLsing the cost fnneiion subjected to the 
■yario'as constraolnts discussed in. section 3*1. The principle of 
opti5:*a;lit3r as s'aggGstod by Bellnan |^8j ,■ is invohed in order to 
deduce a, recurrence relation characterising t!'.-.at fnnetion* ■ A solution 
of recurrence relation is then -sought over decision variable and the , 
o'ptirai decision at ea>ch stare is determined and thus the optimal value 
of the objective function obtained* For the sake of simplicity only 
three disciplines - Mechanical, Electrical siid Civil Engineering have 


been considered- 

Two algorithms are generally used for solving a, dynamic 
prograr.rung . prebien, name ly’ f or^uzrd and' baclcsTard a-lgori thms* ' ' In ■ the' 
forra^rd a.lgoritl:nn'one proceeds ■ from first to last sta.ge-; In^the ■ 
ba, cterard a-lgorithm. one, proceeds from last ' stage,' to the^ first stage , ' 
when 'nethinp us specifie.d' at ' 1-ast. stage - ' It, ,',i S' ' always'" better ; to ,' " 
wo,rk'rrftli'''bacivward;:algorifhm 
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Formlation cf the Lodel : 

Givein the rzininan and riaxinom midber of pernissible se3-ts 
in eneh discipline, dynninic prcrrainrant^ is nsed to deterninG the 
scpnencG cf decisions which optiicises the total expected cost# 

hefininp: : 


rjo 

V 


: the miriber of students likely, to join the discipline 
th 

j for the k period, 

= the nuriber of sea-ts to be provided in the discipline 
t^'' 

j for the k ■'** period 
(decision v?n'‘ia-ble) 

It any pariod k, 7 ^^ is the state of tb.s system for the 
discinliris 3 and is the rsnnber of seats to be provided, for 
that discipline, 'fhers are t”ro distinct possibilities, 

i) the minber of students likely to join the discipline j 

are less than the manber of seats prcTided for the 

J ' G J 

scinline - nndernredictiori Z, 

ii) the niirilcer of students likely to join tne ciiscipliBO' , j 
a,re nore than the nnnber of seats provided for the 
discipline - overprediction ^ 

Let and be the penalties forpander and overpredictioas for • 

" . 1' ^ 2'.„ a'*'. ■ ■ ■ ■ . 

"t ''’"0 di scic^iine !■ ' '.resTce ct:ivo:ly,# ' ■ .Ttiese ; 'penalniGS'' niny , oe or' 

non-linear functions of uj' and . The exact nature of tlieir 
deix-ndonce on V? and sj will depend upon horr severely an 






o o 


total o^rpectel, cost for i!r;dcrprediction Yicnld to obtained hy 

intogratin.3 beto/Gon ioTjer linit and for tlie discipline j. 

Si:rddr'j:ly t'oo total expected cost for cverprediction for the 

discipline j nonid be obtained by integrating between and 

nnper Unit (k!-)'* dnerefore, the' total expected cost of the 

th 

diseioline j for the 1: stage is given by 


ttO 


f ■] 


E. 


,^ij 

1 k 


nr 


k = 1,2, 

3 = 1 ,2, > • # # ,11 


rij /'yrJ _ (f ^ 


nr 


(3.31) 


iTtc rc , 


N :» total mnber of stages taken 
n = to tal nnif 0 c r of . di s c i p i i n e s t '^k e n 

The total ez-oectod cost for discipline j over all stages = 


... 

k=i 

(3.32) 


fc] 


Apnljrir.L£ the -:rinci-',le of opticr-lity, tho tot-T-l c-ctinal cost for tlic 
Dtc>-;e, d (da is expressed ao : 


Si(l4) =, liin 


■ 1 ""' 


/k 




r H -4'!^ 4 i 


CiK.i^ 


(3,33) 


J . s=; ;1.,.::2:,|',^' * 30 ' 
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and the constraints are : 


i) The total nu?:i)Gr of seats to be provided in the institution, at 


any period is hounded. 

Z < ■ h 

i=i 


(3.34) 


wlie TQ , 


= total niiriber of seats sanctioned to the institution for 
the period* 

ii) Tbe msniber of seats in each discipline at any period should, fall 
hetTiFoen its i:iini 1:110 and iraxi nun Tallies ^ 

4 ^ 

where, 

aai are the lower and upper Units on the nudber of 
th 

seats in discipline j for the k period. 

iii) Tho niniaarj and naxinon nunber of peruissibie seats in each 
discipline for four consecutive periods is bounded. 

si Cui’ + uf . +U;’ „ + o 4 

rm'^ k k*fi k^2 . k+3^ nax 


' where , 


S^V and'E'^' ' are the '-ninimn -and ' mximn nun^ of 

' nin : , , .mx ■ , , 


perni'Ssihle, seats' in ddsci.pline'-j. for -four consecutiYe periodst 
' 'The dynaoic' ' prograiim technique yields iiie optinun, " 
number of seats to be provided in, each discipline of tne institution 

,CQnsiderin'g;:the:':'Fabi^eus;:'':Const,:r^hi^a:: l'iBte;^:ahOT 

'.'■is';' c'arrle€v''''0ut;-.;tt's|hg';'the;::'erite^'m 
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'ctrdl s 


the solution procedure used are given below. 


3r4: Solution Procedure ; The steps conprise of : 

STEP I : The set of oatrix equations (3.25) and (3,28) are coupled 
equations and cannot be solved directly. Therefore it is 
necessary to go in for an iterative scheoe. Gauss Si e del 
iterative schene with naxioin pivot strategy seens to be 
fairly good for this purpose. 


STEP II : An initial value for ~p is guessed. Making use of these 

results equation (3,25) is solved to get an estimte for 
(l) (I’i 

Z ' The value of Z^' ' thus obtained is substituted in 
n n 

equation (3.28) and a better estinate for i s obtained. 

The procedure is iterated till the desired convergence is 
satisfied. The values of various variances are assuned to be 
known[27j. 


STEP III: 


STEP IV : 


STEP V : 


STEP VI :: 
STEP VII : 


The stage i s upgraded and the above procedure is repeated, 

A better approxination for vector is obtained. 

The schern is iterated till the desired convergence is 
satisfied, ' 

The whole process is iterated till the two consecutive set 
of values of vector p obtained fr co two consecu tive stages 
natch to the desired accuracy. 

Steps I to V are repeated for each discipline inaependently. 
Values of f, ,fg|g h obtained fron StepV are used 
directly in the set of equations (3 #13) and (3#14) to 
forecast the imoher of students likely to join in future 
periods for each discipline individually# 
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STE? VIII : Koeping in view that these naiiy students ^obtained fron 
Step V’ll'^are now interested to join for future periods, 
en up- or bound on the nunber of seats for each discipline 
is fixed accordingly* 

STEP IZ : The density function for systen state for each discipline at 

each period is calculated using the recursive relation for ' 
variance of systen state » 

STEP X : The functional value in equation (3*33) is obtained using 
composite trapegioidal rule * 

" STEP XI : Finally the optinun niuiber of seats to be provided in each 
discipline and the correspondirg optiiaal cost is obtained 
using iynanic .pro grariiiing* The results of the dynamic 
progrruxiing are given in Chapter IV. The planning 
horizon of five periods is assuned* 

COfIPUTSE 

Corzuter prO'grarxies' hove been' developed 'ana teste ci for ' 
solving all the. three nodels described in this Chapter . A listing , 
of the cooiuiter progranm. is .-given in -.Appendix C, 


(^■LOTBE IV 

isisuLffi mu Discussions 

4.1 lilTDDB'JCTIOU 

O’he optinun miaUer of seats to be provided in the various 
cngincoririi"; disciplines have been obtained using dynamic prograimaiag. 
The optimization was carried out using the criterion of ninimizing 
the cost fi’nction subjected to the various constraints as discussed 
in section 3,3. 

4.2 DC'lSCASm 

4.2.1 Linear Cyclic Forecaster 

The linear cyclic forecaster developed assumes a linear 
groxirtb rato. Three types of seasonality factors are imposed on the 
single demand trend, 'The forecasted demand for the number of students 
who joined the various disciplines i s calculated using the past data. 
Figures 4-1 through 4-3 indicate that for cay academic year the 
nuEiber of stuikjnts .admitted for the first seffiester is always higher 
than the niir.i)er of students admitted for the second semester. 

As the decisions are based on the forecasted values and the 
inaccuracies in forecasting may lead, in sone cases , to unprofitable 
policies. Ee suits indicate that the linear cyclic growth for the 
forecasted demand yieldis, very'-litt^leDdifference; bpt^en^the: actual ; 
and the forecasted values for the number of students who joined 
the various disciplines. 









;'w‘?v^^’?''VV'.' 

Wef&{b= 


s/ftS'' t *t.^L53L-s 




FlG.4-2 ACTUAL AND FORECASTED DEMAND 
- STUDENTS WHO JOINED ELECTRICAL 





41 


As has already been discussed in section 3,2.1, the 
correlation and spectral analysis is used to check the effective- ' 
ness of the iterative procedure enployecl. In this technique of 
cori’eiation and spectral analysis the value of lag nunber (d) 
is varied in the range cf 1 to 8 by following "window closing 
procedure" and tac subsequent values cf correlation coefficient 
and snooth spectral density function (G^) are conputed. The 
values of r^^ and obtained for all the ranges of n are analyzed. 
For significant values of la, the values are tabulated in Tables 
1-A through 1-C. 

Figures 4-4 through 4—6 reveal that for 95^ confidence 
level, correlogran of past denand data does not indicate any signi- 
ficant periodicity. Figures 4-7 through 4-9 show that for 95^ 
confidence level and naximn lag equal to 8, the spikes in the 
povrer spectra indicate the periodicity of 4 for nechanical 
engineering, periodicities of 3 and 2 for electrical engineering 
and periodicity of 3 for Civil Engineering. These values of 
periodicities are again cenfirned with the iterative schene 
(discussed in section 3,2, l). The forecaster paraneters as obtained 
by the iterative schene are tabulated in Table 2. 

The results of the iterati'^e schene show that the total > 
error in the forecasting for the past lO periods is -4, 5 and 1 
unit over the total actual dfesand values of 224 , 248 and 186 for 
nechanical engineering, electrical engineering and civil engineering 
re^ectively, 
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TABia 1-A 

COr-jEUSIClT il-ID SJECTlLtL illAJZSIS FOE l/iSCI-Ifi-IIGAL 
ENGIl-lHEHING DISCI PLII'E 

i'laxinQia lag (m) = £ 


I 

S.lIo.| 
1 

' Eag ] 

i V 1 

L_ _J 

1 Correlation 
\ coefficient 

Smooth spectral | Frequency 

|Density Function | f 

J :■ 

1 

•iji 

1 .0000 

4197.16 

0,0000 

2 

1 

-0.0745 

1725.77 

0.0625 

3 

2 

-G,5002 

14.67 

0.1250 

4 

3 

-0.3832 

142.13 

0.1875 

5 

4 

0.9210 

338# 82 

0.2500 

S 

5 

0,0107 

70.13 

0.3125 

7 

6 

-0.7023 

5.00 

0.3750 

8 

7 

-0.9041 

69.95 

0.4375 

9 

8 

1 ,0000 

135.34 

0.5000 


TABIE 1-B 

C0S;3L/.TI0!I ilD SIEC'rPXL il'lALYSIS POE BLECTEIGAL 
EN(rII®ESENG DISCI PLIIE 

Mazima Lag (m) = 8 


I I Correlatioia | Smooth spectral | Freq^uency 

S»No.v Lag I coefficient | Density fTia.ction| f 


I k I 





1 

0 

l.OOOC 

5223.93 

0,0 

2 

1 

-0.6697 

2225.81 

0,0625 

3 

2 

0.S514 

16,05 

0,1250 

4 

3 

-0.6458 

6.27 

0.1875 

S 

4 

0.7841 

12.30 

0.2500 

6 

5 

-0,7251 

24.50 

. 0.3125 

7 

S 

0,9423 

, 58,34 

0,3750 

8 

7 

-0.9530 

227,55 

0.4375 

9 

8 

1,0000 

, , 100.00 

0,5000 
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TABLE 1-C 

COBr.ELLTIOll SPSCTBAL MAUSIS FOE CIVIL 
ENGIIEBEIIJG DISCIPLIl® 

MaxizaoM Lag (la) = 8 


S,IIo.| 

1 

1 i 

^ i 

1 

1 1 

\ Correlation | 

\ coefficient | 

! j 

Smooth spec- 
tral Density 
Function 

k 

1 

1 

1 

Frequency 

f 

1 

0 

1,0000 ■ 

3923.61 


0,0 . 

2 

1 

-0.5843 

1234,93 


0,0625 

3 

2 

0,24S5 

5.00 


0,1250 

4 

3 

• O.llSl' 

16,38 


0.1875 

5 

4 

0,2814 

79.80 


0.2500 

6 

5 

-0.1351 

206.14 


0,3125 

7 

6 

-0.1617: ; 

23.25 


0,3750 

8 

7 

0,7254 

337,76 


0,4375 

9 

8 

-1,0000 

303.49 


0.5000 . 










FIG- 4-8 SMOOTH F 
NUMBER 0 
.. r DlSClF'LiNE 






51 


?i3IS 2 


?^r. o'rTT:p'?. r":- n Tjrn-niT- o 
3 jX.v h..Jr*i lib: .*,j 


”1^ “/allies for liecbanical j 
Ferai:*eters X Si\iin-3ering discipline] 

f ” ] 

J V&lies for Elec- j 
( trical Engineering] 
[' disciiJline 1 

\ Values for 

1 Civil Enginee*- 
^ rini^ discipline 

A 23.38 

15.90 

33*29 

3 -C.24 

1 #38 

--2.76 

D 5.01 

-13.08 

4.97 

S 4.5 

1.50 

2.50 

p 8.71 

- 8.09 

6.66 

G 2.00 

■ 1.00 

4 *00 

H -3.29 

9.18 

5.03 

U 1.00 

1.50 

1.50 
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4.2.2 Kalrian Filter Forecasting Model 

Vhroe ty',wp of problems - smoothing, filtering and 
proiviftiun as discnssed in section 3.3 can he handled, hy this 
iijodel. ^ocaiiae of its filtering property it filters the past 
states at each succeeding period and therefore the peet states are 
sraoathefl to a great extent. 

Tills is a probabilistic model. The principle of maximam 
likelihood estimate is used for •;rGdicting the optimoa values of 
syston stetos and parameters (discussed in Chapter III). The 
Gauss-Siedol-Iterative Schene with maximm pivot strategy is employed 
for tlio evaluation of the system states and parameters respectively, 
App,endix B gives the flow chart elucidating the Gauss-Siedel proce- 
dure, The procedure converged in 10, 8 and 9 stages for mectonical 
engineering, electrical engineering and civil engineering respectively. 
The optinal values of Kalman-Filter parameters have been tabulated 
in Table 3-A and used for future forecasting, 

Kalroan Filter nethed gives a recurrence relation for the 
variance of the system state at each period. The total error in 
forecristing for the past i0 periods is -3, 2 and 1 unit for 
meelinnical engineering, electrical engineering and civil engineering. 
The forecasted values of the system state for both models for 
planning period have been tabulated in Table 3-B, 

The total error in Kalman Filter model is less than the 
cyclic forecaster and therefore the values obtained by Kalman 
Filter are adopted for pi aiming. Purposes, ®he pr^ 

model is therefore developed on the basis of Kalman Filter. 
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T/iBIE 3-A 


?/iHJ\JOTSPJ3 C-f' 2AIMAN FIIiTER 


I '- .luop '!:i' 'tocIia-|/'E.l'a.es for Elec- | Raines for Civil 
nice.1 BxigineGringltrical Engineering! Engineering disci- 
fliscinl iae ^disein line • | pline 


.8oe 


0«566 




-0.235 


0,035 


0,1T8 


0,317 


0.158 


i 


1.746 


21.418 


2.513 
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TABLE 3-B 


T'e ioieb c? S’rrj^ii'iTs Liiisi;/ to join 
lOL'S 3Y CYCLIC FOEEG/iSTER MD KAUm- 

PILTSR MODEL 



OR'I-FILTHR 


iCKjct.ri cnl 


Civil 


discipline 


Slectri- f Civil 


Bginecriiu Eiigineoring Ccil Engg# f ■ cal Engi— | Engine erii 


discipline discip- 



if neering j 

i di£cipHne 


discipline 


1 

23 

27 

29 

37 

45 

25 

2 

30 

50 

31 

50 

60 

40 

3 

11 

30 

19 

35 

40 

35 

4 

14 

46 

29 

55 

65 

45 

5 

13 

34 

27 

31 

35 

30 
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Ih'.,' iiiizr.tion ;[;rocadurc of dynardc prograoiaing attenpts 


Vv: .t the ciptiimus nn-’ber of serits to be ”rovid« 


io in each 


rlisciplino. The of Kalnnn Filter forecasting nodel is used 

as t!v. i'v ;•'.{■ to d.Wini'ic oroyrarrdng nodcl. 3ven and odd periods 
arc,* v; ci ' i.ir t’K; first and second senoster respectively. 
Tadir:,: the •iri/'ut data for problem ns : 


i) TJ-a rtaximui number of pernissible seats as sanctioned by the 
l.i.T, Kr.’vur a’-ithcrities for ^lechanical, electrical and civil 
C'ligiiiL'C ririf disci pliiies nre 90, 110, 80 respectively# - 

ii) The rini'iu!:! nu;.d)cr of students needed for any discipline to 
exiot is an ach inistrrtive problen which is a function of both 
intan, 3 ible and tangible factor. In this work the lo7/er Unit 
cf 7C, 90, 60 has been assuned for the nechanical engineering, 
electrical engineering and civil engineering respetively. As 
a .’.r.iter of fact ainy values for the lower linit can be assuned, 
Frori the past data, it was felt that the Units stated above 
are m’ite I'or.sonnble. Further as the Sfference between the 
upper arid lower Units on the pernissible seats for the various 
disciplines increases, the conimta-tional tine increases. 

The optiraal policy and thus the optinal cost depends upon 
the penalties inposed. It has been assuBBd that the factor -y 
P-; :■ ( fy « C^/C 2 )ican ■ take^ any valuehfi?©a''0’*l onwards. :/To^: stu^:-- the • 
effect of y on the- optirjal policy optiinization has aeen carried 
out in the range of 0.1 to 10, fable 4 gives results for 
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f ABIE 4 


SFFi3 .:7r 


aiSE OF C^/Cg 3ATI0 


ON THE OPTIIvI/iL COST 


Br.i 


.lti< 




n<K2ring discipline 

S ) = c / G ^ 

X ‘iU 


Vc-lus of Objective function in 
terns of units of cost (optimal 
cost) 


0.1 

O.R 
0,4 
0,6 
0.8 
1 ,S 
1.8 

2.4 

4.5 
10.0 


0,7528 

0.7034 

0.6541 

0.8048 

0,5554 

0.5061 

G.4567 

0,4074 

0.3469, 

0.2312 
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1 r 'if'.', r ' r.t r :: , 

■■■r K’. the rc.'i.ult.'), it is cler.r th.'iit the optimal cost 
■:? or.- r.;-r:.ly as the y-Muc of / increases. This is evident 
s:i'!Cv 1 . ere; tsir-i' iY.i;ueo3 the region of underprodiction, as 

1 v'l . .1 t.i e.. arc i.T*,,uROd for uuderpre diction. This results in 
r,,.:;rr'ir.c!iin,t cUjtcr and closer to 1^. 

T i c ncUr’e, it is observed that r. value- cf -7 equal to 
10 :a;iy iso m ontimal sns for the system turlsr c ansi Oration , The 
o-itiral ralicy I'-t the nunbor of scats to be pr.ovided for various 
disciplines c'Trcarxsnding to > equal to 10 has been tabulated in 
Table 5. 


The optirrtl policy for the nuciber of seats in the first 
oc'aestor of 1973-74 as obtained by dynanic progrcuaiung for necha- 
nical erjftintxsring, electrical engineering .••nd civil engineering 
are 26, 35 and 24 rc*sj:3ectivi>ly. The actual intake however rras 20, 
29 an;’ ■»5 rcspoctively. This discrepency in the results nay be due 
to functional forn of the penalties assumed. The functional fom 
of the 'X'naltioc assumed was linear while in reality it my be 
ncn-li nr. 
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TiSLE 5 


„ , iit(' 'vii 'K,a^ *ji,i .. ,«'..i'„ '^,,<^ 


■5Yiv,nc JTjCv/nffiTG 
^ = ic 


TS'jtinrl I'd icy for fOptinal rolicy for" 
'*' ::'T.nicr.l En.'~i^®e-felectricel Snginee- 
v.': r'isci-rlir.g fring discioliRe 


r~0ptim£i Policy for 
I Civil Engineering 
jj discinline 


1 

15 

18 

15 

?. 

?.3 

35 

24 

3 

IS 

20 

11 

4 

R5 

36 

25 

5 

■«5 e 

r* '*>* 

19 

13 


OgAPTER V 


GONCJiUSTONS AND SCOPE POE PDTCBE WOEK 

CCIICLIBIONS 

In the past, the probleias pertaining to educational 
planning have been tackled using intution and experience. This is 
probably due to .the fact that educational systems are veay complex 
and difficult to handle. With the advent of high speed computer and 
optimization techniques. It is now possible to solve some of these 
problems using modem quantitative techniques. In this thesis an 
attoapt has been made to solve the problems pertaining admissions 
to an educational institution using such techniques. 

The proposed models have been validated using the data 
furnished by the Post-graduate Office of 1. 1. T. Kanpur. The total 
absolute error as found by Linear Cyclic forecaster for Mechanical 
Engineering j Electilcal Engineering and Civil Engineering was 8,C^ 
3,6^ snd 10, lespectively. The error obtained is small and thus 
the proposed forecaster is sufficiently accurate. 

The author would like to point out that this research 
work should be considered as a beginning towards the use of modem 
aanagca-ieBd> tools for educational planning and lot more needs to be 
done to bring in the necessaiy so^istication and reality into 
the proposed models. 


60 a 


SCOPE FOR FUTURE WORK 

There are several avenues for the extension of the 
present work, hsx atteapt nay he nade to relax some of the assump- 
tions nado in developing the model. Because of the ignorance of 
the exact functional form of the penalties the author assume it 
to DO linear. However, it would he interesting to find the 
exact functional form of the penalties. Further the proposed 
model doen not account for the students who leave the institution 
in between the terns and those who overstay. To nake the model 
more realistic, one should study the effect of this factor on 
the decision policies. 
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APTONDIX A 


CYCLIC FOBBC/i^TET^ ; 

In tae linear cyclic forecaster the values of the 
coefficients A, B, D, F and H are confuted using the least squares 
technique# In this technique a function is obtained by interpolating 
a set of given data. Let the function b.e of the form 

JD(t) . f(t) (A-.1) 

where JD{*t) is the forecasted value of the demand* 

The function JD(t) should he such that it minimizes the 

standard error of estiaate. The standard error of estiaate.S. , ia 

t 

defined as : 

N 

^ (Error in forecasting) (N-f), when f (A-2) 

1 

where f = degrees of freedom 

N = the number of available past demand data points ' 

To miniaize S, is equivalent to minimizing E, 

where E = (Error in Forecasting) (A-*3) 

.1 : 

This amounts to minimizing the sun of the squares of the difference 
between the actual value of demand observed at time period t and 
the value t. rat would be predicted from the forecasting function, JI>(t) 
The expiession for E given in equation (A-3) for linear 
cyclic forecaster is : 

N . . „ ■ - . - 2 

B = jj[D(t)-A-Bt“I»(Cos( 'nt/fe))-E(Cos( TTt/s))-H(Cos ( 7TtA)^ 

1 (A-4) 


wnere 
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= the constect which represents the demand at time period 
zero ("y” axis intercept)^ 

I.'i(t;= actual demand for period t, 

3 = slope of "Straight line" 

3, 3,3: <= coefficients which determine the aii^iitnde of first, second 

tiiira. cyclic trends respectively, 

E^G,U = coefficients which determine the periods of first, second 
and taird cyclic trends respectively. 


Tae forecaster pa-rarieters are obtained by minimizing S v/ith 
respect to 1, B, D, E and H. This yields. 


A 




Z [lD(t)-. 


•A-Bt-B(Cos( ^t/E))-P(Cos( ^ t/G))-H(Cos( Rt/U))J=0 

(A-5) 

||= £ t[lD(t)-A-Bt-D(Cos( ?\t/E))-P(Cos( H t/G) )-H(Cos ( T\t/u))J\.V 

(A-6) 

^ Cos( /r t/E) f'ID(t)-A^t-D(Cos( A VE))^(Cos( fi t/G)}- , 

H(Cor( /Xt/j))'] r= a (A-7) 

I'l ^ 

» ^ Cos( 7(t/G)[ lD(t)WL-Bt-^(Cos( rU/l))-F(Cos( 


1 

N 


H(Coo( TTt/u))] = 0 (A-8) 


Ai 


^ *= ^ Gos( IX t/G)| ID(t)-A-Bt-4)(Cos( 1\ t/E))-^(CGs( K t/G))- 


h(cos{ TWdI')] = 0 (A-9) 
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THIS IS A FCRLCAST.R prograb 

PURPQSr - TC G .OfHAT'- FORECASTER PARAHtTERS.* 


PRQGR&MM .. FCR LI NA-iR, CYCLIC D-MANO FORECASTER 
*** DtSCRI PTlOi'., OF Input parameters - 
N = DM=TLiTaL no. of past peri CDS 
IO( J ) = 0=4CTL/,L OEMAoD l:i PERIOD J 

=:«*=!< DESCRIPTION OF VaRIASLES - 

JOU) - FORLCASriR DFMAMC FOR PERIOD I. 

ISRR(I) - ERROR IN FORLCASTING FOR PER! 00 I. 

IirRAB(I) - AOSCLUTw ;RRORIn FORECASTING FOR PERIOD I. 
ItiRTGT - TOTAL .RROR IN FORECASTING. 

JDTOT - total FuRfcCAST DEMAND. 

MINABT ~ MINIMUM ABSOLUTE c.RROR IN FORECASTING. 

NN,MM,LL - JO. OF PERIODS IN THREE CYCLES (DIFFERENT). 
N'N3PT,MMCPT,LL0PT - OPTIBUM NO. OF PERIODS IN EACH CYCLE . 
SIGMA - STANDARD DtVIATICN. 

IDTrjT = f)TOT = TOT AL OF ACTUAL CcBANOS ID(J) 

ROUTINE USED * 

ORCS'^^^ ' 


3^ sjc )}C ^ Jit ^ J^c ;.Je '3^ 5jc 2^ i?c ^ ^ 5{S :3j« 

DIMENSION IC(IOO) ,JD{100)tIERr<(100),A(10,10)f 8(10,10) ,IERAB{iCO) 
READ 200, N 
2Cu FORMAT (13) 

READ l-D' , (ID ( J) , J=„,N ) 

100 F0RMAT(10I6) 

IF(ID( i ) .cC.i )G0 TO 5DG 
PRINT iA}, ,N 

l.i FORMAT (///////iaX,*TOTAL NO. OF PAST DEMAND PERIODS = *,18////) 
MINABT=0 
NM=5 : 

NM1 = NM + 1 ' '■■ 

*** FOLLOWING LOOP CALCULATES TFt OPTIMUM NO. OF PERIODS IN EACH CYCLE 
FOR MINIMUM TOTAL ABSOLUTE ERROR IN FORECASTING. 

DO 2C NN=2,K : ' • 

00 20 MM=2,K 

IF{MM.EQ.NN)GO TO 20 : \ : 

/ • DO 19 LL=2,N 

IFTLL. EQ.MM)GO TO 19 

*** CALL FCRCS SUBROUTINE TO CALCULATE TOTAL ABSOLUTE ERROR IN 
:^D::;poRECASTiNG. 

CALL F0RCS(N,ID, JD,ieRR»IERAB,A,IDTOT,ODTOT, lERTOT , 
iAA,BB,D,E,lA8T0T,H,NN,F,G,MB,LL,U,NM,NMl) 


STfJRc HliMIFlM V/LU,: GF TOTAL ;..RROR IN FORECASTING AND NO. OF 
PE KI 00 S . 

IF (MIN HU. H. .0 

IFCMIN^bT-IALTOT):, , , , . 

1. MINA3T=I ; 3TCT 
NNOPT = iiij 
,MMQPT = ''’M 
LLCPT=LL 

19 CONTINUl 

I’O CONTI HIE 

CALL FORCS SUBRCUTiN TO G-NERATE FORECASTING PARAMETERS FOR 
MIMINU« TOTAL ABSOLUTE CRRCR IN FORECASTING. 

CALL F!jRCS(N,I U, JD,Ii;RR,IERAB,A,IDTOT , JDTOT, lERTOT, AA, 8B , D , c » M IN 
: ABT, Ht TriOPTjF, S , '^MGPT , LLOPT ,U , NM, NMl > 

PRINT 1} :,NMjPT,MMQPT,LLOPT 

112 FORMAT (////3X,') ( .H*)////SXt.*MO. OF PERIODS PER CYCLE (NN) = *,I8 
i////8X,*NC. OF PERIODS PER CYCLE (MM) = =!=, I3////8X,*NQ. OF PERIODS 
2 PER CYCLE (LL) = I6////PX, CO( IH*) //// ) 

PRINT 1u 2,AA,8B,D, r,F,G,H,U 

12 FORMAT (//SX,4 THIS IS THE DEMAND FORECASTER *«**////5X,63 

IHD (T ) = A + B*T + D*COS { PI*T/ E ) + F =>CGS ( PI*T/G ) + H*COS ( P I *T/ U ) +VAR// 
2//8X»’®'A = j, ij • 3 j X t » 

3fcl8»3»5Xt*0 = *jt. lo.d/dX»’5'b = *fc^8.3»5Xf*F=*»El 
48.8»5X»*G = c '.8 . 6» 5X» = * > £i8 « 8/ 8X» =*0= =S‘»£i8.8//) 

PRINT i:3,(I,ia(l) ,JO(I ),IEPR(I),IERAa(I),I=l,N) 

103 FORMAT ( //3X,4 :3H***THIS I S DEMAND FOR PRESENT *N» PERIODS ****#/// 
i/3X,*PtHiOD ACTUAL CALCULATEO ERROR IN ASSGLUTE ERR0R*/15X, =i'DEMAND 
2 DEMAND FORECAST IN FORtCt^ST*//// (4X,5I8// ) ) 

PRINT i.v'-i,ICTOT, JDT.jTiI£RTOT,mINABT 

104 FORMAT ( //3X , ^TOTAL ACTUAL DEMAND = * ,l 8// //8X , «TOTAL FORECASTED 

IDEMAMD = =4=, I9////8X,«T0TAL ERROR IN FORECASTING = 18////8X t*TOTA 

2L ABSOLUTE ERROR IN FORECASTING = *,187/) 

THIS CALCULATES THE MINIMUM VALUE OF STANDARD DEVIATION. 

^;{e *!{! 4:^! # :J« :{e # aH 4 4 * 54* 4 ❖ 4 * :5e 


lERTOT=IERR(i) 

16R2 = IERR(1)*=4'2 
iER2T=IER2 
DO. 60 ; i=.2»N ^ 

IER2 = IERR{I )*=4'2 
ieR2T=IER2T+ieR2 
ISRTOT=IERTCT+IERR(I ) 
ISUM2=IERT0T**2 


PSlG=FLCu^T ( (l^ra.'iT-ISUM^)/ ) ) 

SIGMi=SQirr{PSlG) 

IF(I.L5,5) GO TO 60 
QS I G =F in;' T ( I -;R^ V ( I-'MH) ) 

SIGMAx = S0!U(0SlG) 

SIGMA^i 01:0 . { S I G‘* . , G I 3MA 1 ) 

60 COiOTIOUl 

PRINT ; 2 CfSIGO \ 

120 FORMAT {3:<,>!-AK = :.a.o) 

PUNCH ^2.. ,AA,iO,D,_,F,G,HtU,SiGMA 
GO TO 20. 

121 format C’XtFio.o) 

5 2 COiMTIHUL 

STOP 

END 



BFTC 


SURR OUT I ,-L- FOR C S { U , I D , J 3 , 1 £ RR , l ftR ab, A , I DTOT, J DTOT , I ER TOT , A A , BB , 
» c » I ^ BT;jT I h >i i ^ j F , G, WM, LL ,U ,M!«! ,NM1 ) 


^ it : 4 < 4 ; ^ 


ii -J|c :{c 3{c Jjc :{; 3j: 5§! 3j: :!{c sjt ^ ajc ^5{c3{c3^5^::5c::§e':^!J^c:^ 


THIS 1$ SUuR'JUTIMS FORCS 

*«* PURPOSE - TC ic MINIf'UM ABSOLUTE ERROR IN FORECASTING. 


5js ^ 3*< if /f it if: 




D I ME NS I GL i C ( _ 0 'I ) , JD ( „.a*} ) , I ERR ( 100 ) , A ( 1 0 , 1 G ) , B { 10, 10 ) , I ERAB ( 1.00 ) 
cNN=NN 
E = ENN/*;. 

G=FLOAT(MM) /:. 

U=FLOAT{LL ) /£. 

IDTOT= 

DO 1 J = .l ,N 

1 IOTOT=IDTCT+ID{ J) 

DTOT = I OTi'T 
DN=N 

DAVG=DTaT/CK 
DO 7 L=i,NMi 
DO 7 M=i,NMl 
7 A(L,M) = . 

REDEFINE D AND T BY AK AND F TO DECREASE NO, OF CALCULATIONS 

DAVG= AVERAGt ACTUAL DEMAND IN N PERIODS 

AK=D-DAVG 

T=J= PERIOD 

H=T-N/£ 


if if ififififif':iiififif:if:if:^ififififififiiififififififif^:^if:^ififififif;^ifififi 


REDEFINE D AND T BY AK AND F TO DECREASE NO. OF CALCULATIONS 

AK=D-DAVG 

T=J=PERIOD 

A(3,l)= TOTAL H*AK 

A{2,4}=A{4,2)=0=T0TAL COS 

A(2,3)= A(3,2)= TOTAL H 

A(2,2)=N 

A{2,1)'= TOTAL AK 

H=T-N/2 

A'{3,3)=A(3r3) = T0TAL H CUBE 

,AI3,4)=A(4,3) = A{4,4)=N/I ; , ' ^ 

A(3,l)=T0TAL OF H SQARE4AK » 

A(3,4)=T0TAL H SyUARE’FCQS=A (4,3) 

A(4rl|:=TaTAL OF K4,CQ:S:- ■ i: ; a: 



Di.) 1 J = * i''i 

r)=iL)(j) 

AK=D-l'';: VG 
T = J 

Cu = :us { { ...£.*T) /( 7 .--i=r ) ) 

C 02 = cns( ( 2 :'.#T )/( 7 .^<G) ) 

CDi>=C:JS( (....=:•[) /(7.*U) ) 

A ( 2 : , 2 ) = A t 2 , 1 ) + / K 
A { 2 » -j ) = w { 2 I 2 ) + 1 1 

A ( J » j> ) = / ( !• f 3 ) + H ■'!' H 
ACi, a,)=A{3,l) + H'-l‘AK 
A { «■ f .!. ) = A ( j 1 ) + H * H * A K 
A { A » n } = A { A f A ) + H H * C U 
A( '.,4) =A(^,A) + CO 
A(3,4} =A(3,4) + H--!‘CG 
A{ 4 » 4 -} = /i {4 f A ) + C i C u 
A(2, 5)=A(^,5)+cn2 
A(?, 5) =A{3,5) + H’i‘CCr 
A(4,5J=A{4f5) + CD*CC2 
A{5, i)=A(5,l) + AK=!'CQ2 
A(5, 5) =A(5,5)+C02*C04 
A( 2,6) =A(2,6)+C03 
A{3,6)=A(3,6)+H*C03 
A ( 4f 6 ) — A ( 4 » 6 ) + C Oj) 

A( 5, 6) =A(5,£) + CQ2»^Ci:i3 
A ( 6, 6) =A ( 6 » 6 ) + C03^C0 * 

A{6,1) =A(6,l) + AK*Cij3 

2 A(4, iy=A«4,I)+AK*C0 
A(2, 2)=UN 

DO 4 K=2,N^' 

I = K 

3 1=1+1 

A{ I,K) =A{K, I ) 

IF(I-NHl)3i4,4 

4 COrjTIMUE 

DETERMINANT OF ORDER NM BY PIVOTAL GQNDENSATION METHOD 

44 44 44:44i444 4 44444^4:*=ii‘^***^**^*^*^****^^*******=i<*^*****^=^***^***^^***’i‘***^ 
DO 35 II=1,NH1 : 

• DO 6 J=1,NM 
:DO 6 : K=i:, NM G ■ 

Jl=J+l 
KI-KE - A 

'IF ( K.GE.n )Kl=Kl+i a: 

6 ^ 8( JrRT-AiJlYtt 

•:e: 1^42: vEy;:::A;:;^ -^VE •: ^a.a : ; : ^ a; a ; r -vA^^'-EE^r 


OcTE = ; . 

M = 0 

3 DO 1 1 = K,;-jN 

1 F ( A BS L } L n . L- . . E-25) GD TO 8 

GO TO .2 
5 M=M+ 1 

DO 9 KH=L,NN 
Li'4=L+''^ 

B ( KM,L )='.: (KW,L -O 
B(Kyi,L''-i)=i3(KM,L) 

9 DETE RM=-Di.TERM 
GO TO 

RATI D=n {1 ,L)/B(L,L) 

Du 1 J=K,DN 

i- 8 ( I , J) =b{ I , J )-B (L t J) *RATIO 
IF (K-MP) ..5,2 ,2' 

15 L = K 
K = K.+ I 
GO TO 3 

2u DO 25 1=1, RN 
25 DETtKM=L)c;TlRM«8(L,L) 

GO TO (3 ,32,32,34,30,37) ,I I 
3e Bl=D£TiRM 

32 Dl=D5TtRM 

33 £i=D£TwR.M 

34 Fi^DETERM 

36 Ol=DET£RM 

37 F2=DtTERM 

35 CONTINUc 
F=DM/2. 

AA = {DX+::**F)/B.+DAVG 
6 B =- E i / B i 
D=+F1/BI 
F=-*Oi/Bi 

H=+F2/51 . 

JDT0T=0 
IERT0T=0 
-IABTOT=i' , 

DO 4a I=i,K 
AI=I 

JDd ) = AA+8B*AI +D=^=C0S( (22.*AI ) / (7 .*£ ) ) +.5+F*C0S( {22. *AI ) A( 7 
I.^G) )4H*C0S((22.4AI}/(7.*U) ) 
lERRd )=JO(I )-IO(I ) 
lERABt I) -IABS( lERR^ 
dABTOT = IABTOT+iFRABrir 
JDTOT=JDTOT+JD{ I ) 

40 IERTOT = IERTCT+IERR(I ) 

RETURN^ ^ 4 ;;;:44 

END 2 - - ' ' E 


;bftc 

4: '-Js -4^ *-5« ^ s;« ijc 5{c 

FORcCASTSR FOR FUTURE H PERIODS 

* * .'It s}: Jjc Jjs :{t ^ ^ ^ ^ J, ^ ^ ^ ^ ^ Jj, ^ ^ 3j. ^ ^ ^ J }5 ^ 


DinENSIC i JC{i ) 

ITR = 

READ ;t.0S,8IM, BI OCR , FLAST 
: 3 FORMAT {3F5. 3 ) 

RcAD 1 4 jSE£i) 
lO-i* FQRHmKF.. 0.8) 

READ 1''>S,N 
li..5 FORMAT (IS) 

1*/ ITR=ITR + ,. 

READ 101 TAtCfE , F , G ,H , U, SI GMA 
lOi, FORMAT (SFltj.d) 

PRINT i ?. , A , D , 1- , F , G , i l , U , S I G ^ A 

I ?2 FORMAT (//3X,*D-MAND PAR AMET ERS ( INPUT ) =i'/8X , 25 { 1H- ) ////8X, 

a 63H JD { T ) =A + B*T + D=i'COS ( P I *T/ E) +F*COS ( PI’i'T/G )+H*COS { PI *T/U ) 
Z ////8X,’!'A = *,F;.6.6,5X,’!-D = =S=,Fl6.6,5X,*e = «,F16.6, 
35X,*F = *,Fl5.u//8X,*G = *, Flc.6,5X,=4=H = *,F16.6,5Xt 
A-*U = *,F,i6.6,5X,*SIGMA = *iFio.6//) 

PRINT iOi tSEhO 
8=BINT 

II CONTINUE 
IFIB.GT.BLAST) GO TO 99 
Nl=N+l 

ITOTD= 

ITOTA= . 

PRINT i 6,3 
DO 79 I=1,M 
T = I"-i 

MEAM=A+B«T+C*CaS { (22 .*T)/ (7 .*£ ) )+F*COS { ( 22 .*T ) / (7.*G) 3 
1+H*C0S ( (22.*T) /(7.=!=U3 3 + 0.5 
ALL1=FLQAT{MEAN3-3.X‘SIGMA . 

AULi=FL0AT(M*-AN)+3.*SIGHA 

CALL JAYA(MEAM,ALLl, AULI,SIGfAA,SEED,DEVT 3 

JD{I 3=OEVT+C.5 

IDIFF1=JD{ I 3-MSAN 

IABSI = IA6S(IOIFF13 ^ 

ITOTD=ITOTD+IOIFFI 

ITOTA=ITOTA+IABSI 

PRINT 107, 1 ,M£aN, JD{ I 3 , IDIFFl, lABSl, SEED 
79 CONTINUE 

PRINT i08,n0TD,IT0TA 


+ VAR 


B = B+BI 'iCR 
GO TO 

99 CONTI^JU!:; 


2 ji FUi<M«T(//u 
1,6 F0KM.'\T(BK, 
lACTUAL AC 
IDIFFCR a.,|Cl 


A»-Iij!TlAL SEcD = *,FI5.lO//) 

^VALiJ’:. ^ CF GRADI FOT B = *, F9 . 5/ /6X, *PERI OD MEAN 
iCAL ASSCiLUTE VALUE' OF*/ I4X » =!=D EMANO DSMAfsID 
CIFF;,R-:vCL SEED*) 


,07 FORMAT { 7 a, 14,3 
« FORMAT {E.<,*TOTAL 
916) 


6!3TiO,3X, .iO) 
D:r.V3ATE =*,I6»8X. 


I6,8X,*T0TAL ABSOLUTE DLVIATi 


IFdTR.L-U. j),STQP 


GO TO lO 





I3FTC 


bU'JRiJU i I "Ic JAY 1 ( LM jTLMI i'jj-TL NAX, STD, XI , DI ST ) 
XA = WNr)Y. {XI ) 

CALL SNOYKX:, ) 

4 RA=RXDY.;(X) 

IF (RA. L'l ... '3- 5)00 Tn * 

RB=RXaY:: {X) 

IF{R3.L:...’0.oi-,..5)GJ JO 4 

^={-2. , «ALGG(RA) ) . 5*C0S ( 6. '33*RB) 

VAR= V^STD 
DIST=Vi,R + FLC'AT(LY) 

XI = R 3 

IF{DIST-TLFlO)6,7,a 

6 DIST=TLMIX 

7 RETURX 

8 IF(D1ST-ILM/)X)7,7,9 

9 DIST=TLMAX 
RETURN 
END 



IBFTC 


4: ;1< j;< 4: * :{t :;s Jj; v 5jt ^ ^ ^ ^ ^ 5j. ^ ^ ^ ^ ^ ^ Jj, ^ ^ ^ ^ Jjj Jj. ^ 

P:^.OGRAH FOR CORRcLATIGN AND SPECTRAL ANALYSIS 

j,? ,;*< ;*5 5»s: 4c jjc ^ 5^ 3j ^c ■;{« sjc jjc 5{c sjc 4c-5*c ^c 4c 4: 4c 4c 4c ^ 4c '5^:4c :4c 

❖ ❖❖ THIS IS PR T A PROGRA ^❖4’4c^4«^2j:4c.4«4;4c5^4i4c 4c4c4:4c 4:4c5}:4t4c-4c4c4c5^3!£4c:^4c4c4}c:^4c4c4c4c 

PURPGStf - TG TEST CATA FOR CORRELATION AND SPECTRAL ANALYSIS. 

*** DESCRIPTION OF INPUT PARAMETERS - 
X(Lfl) = past demand FI STORY 
L= VARIABLE NO. 

1= SAMPLE NO. 

NSAMP - NO. OF SAMPLES. 

FOR AUTOCCRRELATIiN NVAR = S. 

MLAG - MAXIMUM LAG NUMBER (BETWEEN NSAMP/20 TO NSAMP/2). 

*** DESCRIPTION OF VARIABLE - 

TPd) - SUM OF VARIABLE JX TILL PERIOD I. 

FPII) - SUM OF VARIABLE JX TILL PERIOD I+IP. 

SPd) - SUM OF SQUARE OF JX TILL PERIOD I. 

GP(I) - SUM OF SQUARE OF JX TILL PERIOD I+IP. 

CP(I) - SUM OF PRODUCT OF JX TILL PERIOD LAND I+IP. 

RP(I) - SERIAL CORRELATICN COEFFICIENT BETWEEN VARIABLE JXdJ 
AND JXd + !P) . 

WHERE IP=P=LAG NO. {0,1,2, ...fMLAG). 

WPin - AUTOeOVARIANCE FLNCTION OF ORDER I- 
FLPd) - RAW eSTIHATE OF OF A TRUE POWER SPECTRAL DENSITY 
FUNCTION OF HARPOMC NUMBER I. 

UPd) - SMOOTH gSTIMATE CF A TRUE POWER SPECTRAL DENSITY 

AVEG - MEAN VALUE OF RP(I). 

SIGMA - STANDARD OEVIATICN OF RPd). 

FREQ - FREQUENCY. 

CORNU - NUMBER OF DEGREES OF FREEDOM FOR SPECTRAL CALCULATION. 
GFBAR “ MEAN VALUE OF SMGOTH SPECTRAL DENSITY FUNCTION. 

DIMENSION FMT(l"},X(i,lTC) ,TP(100),SP{100),FP{100),GPd«D), 



( J. G^i ) , WP ( 1 Ji) ) , RP ( 

t ,TiTLu 

RhA 3 it nS '4 V p , i'j 7 « j-^ 

;• FUI<M 4T{;, i5) 

IPRQd=. 

838 CONTINlir 
PRIMT 7 

PRINT 50 .,TITL- 
PRIMT 7' 

6.4 FORMAT C Ao ) 

7;.*. FORMAT ( LX, ) ) / ) 

ifcmvar) .’gc,ll-,:' 


) »UP( 1 ),FLP(1'.G) ,TITL£{ 1j) , JX(i. 


2 •: CtJNTrjiJ;, 

READ4, ((/;(L,I) , I=,’, ,MSaMP) ,L = -^ ,NVAR) 
4 FORMAT (;.rF6.,i) 

■N=NSAMP+:. 


AVE3 =- .. . / ( F L QA T ( i-lSA MP ) - ■ , ) 

SIGMA=SORT ( FLOAKMSA'^P)-:-:, ) /FLOAT(NSAMP-f ) 

TWOS IG = , .96>i-SIGMA 
CURUP =AVcG+TWuSIG 
CQRLOW=AV£G-TWaSIG 

PRINT GO, ,AVL-G,Sn7MA,TWaSiG,CuRUP,C0RL0W 
101 format (//2X,’!'PARAM£TERS for C0RREL0GRAM*//5X,*MEAN =*,Fi3,6,=«' 
iD. DE;V. =*,F1,’.6,* TWO SIGPA =*,F10.6,* UPPER LIMIT =*fF10.6 
Z* LOWER LIMIT =*,F; .6//) 

15 CONTINUE 
READ 3, ML A3 
3 FORMATE 15) 

IF {MLAG.LE.O) GO TO -OO 
DO 300 IX=l,NVAR 
DO 300 IY=1,NVAR 
TP(i )=-. 

SP{i)==, 


FPd ) = 0. 

GPE1) = . . 

DO I 6 l = i,,NSAMP 

sp(i ) = sp(i )+x{ IX, I 

TP (I ) = TP(1)+X( IX, I ) 

GP{1 ) = GP( 1 )4X( 1Y,I )«=i»2 
FPd ) = FP 11 )+X( lY, I ) 

16 CONTINUE 
Mi=MLAG+l 
DO 20 1=2, Ml 
G ,J=I-1 ■ 

K=NSAMP-l+2 

TP(I ) = TP1 J)-XdX,J) 

spd )=sp(i-iy-xdx,j)**2 

FP(I ) = FP1 J)-X1 IY,K ) 

GP(I ) = GP{I-1)-;UIY,K)**«'2 


ST 



2 C 1 J * i T I \ Sj : 

£3 '^iLA3 >’'LA‘,y + l 
30 3\ I=. . 

isiMr'jp=':S'.r-?-i + 

CP{I)= . 

u'U 26 J = , jlvf'I’IP 

K.-=J+I-. 

C P ( I ) =C P ( i ) + ( A ( I X , Ki ) ( I Y , J ) ) 

WPCl ) = C?(I )/( FLOAT (NMINiP } ) 

RNUM=F Lij AT ( N'"*! I iP ) *CP ( i ) - FP ( I ) =1<TP ( I ) 
r<0i_0 j - SlKT ( ( FLiuT { l'P*PJP ) *GP { I ) )-FP( I ) ^^2 ) 

RDcM2= SWHT ( ( Fl OAT { ‘'■iPi INP ) sl^SP ( 1) )-TP ( I ) ) 

3' RP ( I ) = S<i JUP/ (RD-OL^ROLN'J ) 

IF(IX-IY) 2C2,30b,ZO-'' 

202 PRINT :.50,IX,1Y 

i5< FORMAT *CR0S3 CLRRLLATI OK A=*,I4f* Y= * , 1 4// 5X t =i'P*» 7X» *SUM X’X, 
*^4Xt*SUH Y*,;4X,*SUM X S 0=!= , il X t *SUM Y SQ*,11X,*C PR0D*ti3X,*R COEF 
2*/) 

GO TO ^04 
2. .3 PRINT 99,1X,1Y 

99 FORMAT ( / / 2 X , «4UTOCORRELAT I CN X=*fI3,* Y=*, I3//5Xt*P*»7X,*SUM X*, 

ii4X»*SUM Y4,*4Xj*SUM X SQ*,n.X,*SUM Y SQ4,11X,*C PR0D4tl3X,*R COEF 
2*/) 

204 CONTINUt: 

DO 5' I = L,MLAG:;, 

IP = I-1 

51 PRINT I' »IP,TP(I) ,FP{n,SP(I ),GP(I),CP{I),RP{I) 

1 . Al . FO RM A T { 1 H , I 5 , 5 ::1 9 . 8 , FI 6 . 8 ) 

HLAG2=MLAG-1 
33 DO 39 I=*,MLAG„ 

FLP{ 1) =• . 

00 3 6 J = ::tMLAG 

36 FLP{ n=FLPn ) + 2. i*WP ( J ) 4C0S B .1415927*FL0AT( I I-l ) «( J-1 ) ) /FLOAT! 
IMLAGJJ 

39 FLP( I)=FLP{I) + WP(i) + kP(MLAGI)4C0S(3.141614FL0AT(I-l)) 

UP(1)=' .464FLP{2)+.544FLP!1) 

0P(MLAG1 )=0.465!‘FLPIMLAG)+.544FLP{MLAG1) 

00 43 I =2 »MLAG 

43 UP(I)=< .234fLP{I-l)+,544FLP(I)+.234FLP(I+l) 

PRINT 98, IX,IY 

98 FORMAT ( //2X j^SPECTRAL ANALVSiS 0F=«',I3,* VS*, 13//* P*,7X,*C0V 

lAR*, 14X,4RAW SPECTRA*, 5X,4SMC0TH SPECTRA*, 5X,*FREQUENCY*//) 
TOTUP=.. 

4' "DO 75 I = i:,MLAGi- " 

"IP = I--1 

FR£Q=FLOAT( IP)/(2.4FL0AT{MLAG)) 

TOTUP=TOTUP'^UP (I) 

75 PRINT lOL ,IP,WP(I) ,FLP( n,UF(i),FREQ 

CORNU= {2.«FLQA*{NSlHP)/FLQAT{:’^LAO0>-2.:/3.- ' 

gfbar=totlp/flu4T(’'^lag: ) 



,CC^ MjGFr^R 

lo: FOR*■1^ T {//:x,*rA?i. ••.: t".;^s for spectral 

i =*, F:,.6 .2// ) 

3 ;..' CIjNTI\U2 
Gi] T J .. i 
400 L UK'T 1 ''J S I ■ 

I PRO?' = i p'.‘,,jtf + ., 

IFd P^!.;.-..LT.' )GO TG cGO 

STOP 

END 


ANALYSIS*//5X,4=NU = 


*,F 16 . 4 ,>!= 



; :i|: ^ >!« ? 


.j:«sSs4c:s:«55(4«:^< :$5*s:;{; :^:5:‘3jc.i}::{:5Sc35e34<.^3ic35<5:: :^ic 


=>RaGiU ’;■ FGR KAlf^AM FlLTcR ANO DYNAMIC PROGRAMING 

>it * ^>J: ^!s!« **;}!;}: :i< 4 i}!ii: 4;** ;{« #4:* ❖=5'’!'=!'**^** 

vain prog^mmya 


>i;{!ij;sjc:{£5;c;}!!(tjjt!{!:jc;i[>{!sjc!!s#3it:{cs!«!{!:ic^^#^>!'*^!5'=i‘^=i' 

N FILTERING TECHNIQUE 

T£ 


NY**Y**=!'-! 0 . of STAG;;S 
p ( I) *^c:^:=!^PARAr■1£TGRS OF KALMA 
1,(1) *«**SALARY GRADI GNT 
Y { I ) « 0 Li S S R V Ai T I 0 M S 

SI GV ( ! ) *VAR I GNCG OF THE STA 
SI GW***-‘YAR lAMCH UF MuISE IN GBSERVATION 

ERROR BSTW'EN TWO CONSECUTIVE VALUES OF 
IN A STAGS 

cRROR BtTW:::5N TWO CONSECUTIVE VALUES OF 


pjH4:4«**#4:FRACTI GOAL 
PAr<A'!>T'=RS 
HK*«****FRACTi GOAL 


PARAMETERS IN CONSE 


P( I ) «#«#(jbtainl-d for two ST 

X M E C H* >!'=>= P^, A X I Pll) H LIKELIHOOD 
X£LECT^**MAX1MUM LIKELIHOOD 
XCIVILY=4'XAXIMU'1 LIKLIHQOD £ 

5I(5my*«*variance of state 0 

SIGE****VARIAMC'- OF STATE 0 
SIGC****VARIANC: OF STATE 0 
Cl AND C'-' PEMALTISS FOR UND 
Q( j ) ij<*#^HAVE THN USUAL MEAN 
B( It J) =«'*HAVE THE USUAL MEAN 
Odf J) ’!«=FHAVE TH^: USUAL MEAN 
R( I ) ****hAVE the USUAL MEAN 
sue ROUTINE MITA EV 
«#4,«*:«t4**SUBRrjUTINE RAJA EVA 
:{e!5t;^!#>}!s(!’>:*$UeRUUT! fiS INVERT I 
#«*»*«;5:ji=SU8BOUTINE OYNAH US 
TO CPTIPIALITY PREDI 
VARIOUS DISCIPLINES 
:{c#55t#5isj}:**SURROUTINeS FUNliFU 
FUNCTiaNAL VALUES 
*#=^!4 c#««*SIJBROUTINE AVtRG CA 
** * S U B R OU T I N E C OcF F E V 
FUTURE STATES 


CUTIVE STAGES 
AG^S 

ESTIMATE OF STATE OF MECH. DISCIPLINE 
ESTIMATE CF STATE OF ELECT. DISCIPLIN 
STIMATE OF STATE OF CIVIL DISCIPLINE 
F ■''.eCH. DISCIPLINE 
F ELECT. DISCIPLINE 
F CIVIL DISCIPLINE 
ER AND OVER PREDICTION 
ING when reference TO EQN( 

ING WHEN REFERENCE TO eQN( 

ING WHEN REFERENCE TO EQN( 

ING WHEN REFERENCE TO EQN( 

ALUATES THE MATRICES B AND Q 

LU ATES MATRICES D AND R 

NVtrRTS A GIVEN MATRIX 

cS DYNAMIC PROGRAMING TECHNIQUE 

CT THE NO. OF SEATS TO BE PROVIDED IN 


OF 

OF 

OF 

OF 


TEXT 

TEXT 

TEXT 

TEXT 


N2,FUN3iFNl,FN2iFN3 EVALUTE 

LCULATES THE FUTURE STATES 
ALUTES THE VARIANCE FOR THE 

:f jj£ :«! «:«C9{= =5= ********* ********* 

si# ;^ ******* * * ******************** 



IB FTC 


CALL PHJi{32 . ) 

CO^vMU».,/T J/U,Y,CIGV,SIGW 
CQvm0N/IITK/C1,C. 2 
C □ U \1 / 1 i T / 1 SC A T , I r-'' I TA , iMNS 

Ca-'^M0\/,U^.;.'.1/X •, .CH,SIGM,X£LECT,SIGE,XCIVIL »SIGC 
01 YEi-iSIC;,' ISi.AK?) ,IYITA{2,6) ,NN3(6) 

Di i^EnSl La iJ( i.O ) j Y( XIO »SIGV( 2) ,F{40) ,Q(40 ) , B(40»40) , 
..FF(4),i)(-. ) 

DIME^4S10, R{-i- ) ,FFF( :,2,4),X (4., ) 

Oi YE MS ILL! XYEChI.O ), SIGM(.lO ),;;5LECT{iO) ,SIGE{ 10 ), XCIVILliO) 
DIMt MSI 1.1:1 SIGCC 0 ) 

DIMf'jSiGM SIGC, ),KX(i>) 

M=4 

C,^=0..^ 

C2=I.i 

it=l 

3. 2 rf=iT+: 

N=iO 

IFd T. i-0. 

GO TO (3 .3,3^ 4,3 5), IT 
3m3 PRIMT 3'.,6 
GO TO I.O.. 

304 PRIMT 307 
GO TO 3: * 

3i:5 PRIMT 3' S 

306 FORMAT (/50X,=S‘F0R£C'f5TING OF M.5CH. ENGG.DEPTT. */50X, 37 ( iH* ) / ) 

307 FORMAT ( /50X,=5'FORcC,:. STING OF ELECT. ENGG. DEP TT. */5tX , 37 ( IH* )/ ) 
3 .3 F0RMAT(/5''X,4F0RECSSTING OF Cl VI L . ENGG.DEPTT. */50X,37 ( iH* ) / ) 
3i. l READ23 , ( U { I ) , I =1 , N ) , ( Y 1 1 > , 1= . ,N ) 

200 FORMAT! 8F10. 5) 

ReAD£:*.l,SIG^,SIGV(I) ,SIGV(2 ) 

2Ci F0RMAT(3Fiv.5) 

READ212, (F(I ), 1=1,4) 

202 FORMAT (4F10. 5) 

PRIMT25^, (UM) ,I = 1,N) 

PRIMT25 : , (Yd ) ,I = iiM) 

PRINT251,SIGW, SIGVd ) ,SIGV{ 2) 

PRINT252, (Fd) ,I = if4) 

I5C: FORMAT(lX,i.2F23. .5) 

251 FORMAT(i:,.-X,3FI5.5) 

252 FORMAT (10X,4Fi5. 5 ) 

■N=i: 

. " ;^MN=M+l 3 . 

:00: . 3;;: 3 3 v M 

^ QtMm»0 , ■ : 3 : 



DO 1. J = : , A 

1 ( I , J) . c 

',)ij ?: 

D Cj t J “ f ' i 

R { 1 ) = . 

2 

no 3 ! = , , ^ 

3 FF( I ) = F { 1 ) 

6 o C A L L '^1 X T • « ( fv K } L's j y f F ) 
yr4=2#N.. 

CALL ! ';V..nT( '1''i,3,C,X) 

I 

CALL raja (N,!,ij,R,X) 

IF-(O.L'. )3C TO :iO 

MM=4 

CALL I:..V :FT{W‘j|,o,r,p) 

PRINiTl ' iO» (F ( I ) , !=■ ,4) 

300 FC'RMAK :nx*ST ’>Gu^.,OX*VALUE OF FF*/20XtI4, 10X,4( F15.8, 5X ) ) 
DU 60 I = 

HH=( F( i )-FF( i ) )/FF(i) 

IF<A8S(HH3 .GT. . :-)G0 TO 63 
60 CONTINU_; 

DU 65 1 = : 

FF{n=F(n 
65 FFF(N,n=F(I) 

PRIMT iOQ, {X{I ) , I = ID 

lT . format ( 5 X,=» value OF X={=/S (1 X, Fi5.8,5{ 5X ,F15.8) /) ) 

N=N+1 

NN=N4l 

IF{N.GT.iO)GO TO 5 
IF(M.eU.2)GC TO 66 
D067 I = 1 t M 

HK = { FFFdM-: ,n-FFF {N-2t I ) )/ FFFdM-ltl ) 
IFCABSIHKJ.GT.O.ODGO TO 68 
67 COMTI'MUe 
5 C0NTD4UE 
Nl=M-i 

PRIMTlGOOf Nl* { FFFCNIVD rl = l »4 ) 

CALL AVERAG(X,XX,Nl,FFF,SIG) 

GO TOI^P' 3,4C4,4 5) »IT 

403 00 406 1=1,5 
SIGMd )=SI6{ I) 

406 XMfcCH( I )=XX(n 
GO TO 4 La 

404 DO 407 1=1,5 ^ 

^ 

4hJ 7 xeLECT(n=XX{I ) 

AAV'^: GO do 411 lA: 0:1 n :0 : i;;,! -Oo I-: ' 1 ■ i • ' : 

405 -00 408- 5 ^ ^ ^i: 01 

SIGC (I j = SIG{I) 




IB FTC 

SUBflGLirf, i’. fA{NN,B,Q,F) 

C U / T J / U f Y , s I G V , S I GW 

I) I '"Y t'!S 1 1,,;M U ( ) , Y ( 4 . ) , S I GV ( 2 ) , S { 40 , 40 ) , Q { 40 ) , F ( 4?^ ) 

11 = 0 * 1 ;'; 

H*=l .O/ilGy ( i ) 

H = 0 • / S 1 15 \/ ( I ) * * Y 
H3=i . .VSiGw*^=0 
H4=F ( ) ** 0 

H3 =F { } **0 
H6=F ( ■ )*F{::;) 

H7=F(4)**2 

00 : i=:,ii 
Q<J i J = t I i 

1 a( 1, J) = . 

B( 1, 2) =-H6*F::: 

0(2, i) =ri{:.,2) 

B( 2, /.) =-H3*H :-H7*H5-Hi 

O(^) =F{0 )*F(3)*U(i)*H2 

Q { 2 ) =-F ( f ) *Y (I )*HO+F(2)*F{OV)*U{l}*H2 

B{j,,B)=0.0 

3(1,4) =F(1 )*H2 

B { 2 , 3)=H .:, 

6(2, 4) =F( J)*H2 

1 F (ON. i- 0.2) GO TO iO 
DO 2 . K = :i,N 

J = 2*K-;: 

JJ=J-2 
JM=J+I 
B{J,JJ)= . 

B(J,JJ+1)=H1 
B{ J, JJ + Z)=-F’;,-H4*H2 
B{ J, JJ+3 )=-F6*l-i2 
B( J, JJ+4)= • 

B( J, JJ + 5)=F(.l)*H2 
B( JM,JJ)=F(1)*H2 
8{ JW|,JJ+:U=F{2 )*H2 
B( JM,JJ+2)=-H6*H2 
8( + j )=-(H5 + i. )*H2-H7*F3-Hl 
B( JM,JJ+4)=Hi 
6{ JM,JJ+5)=F(2)*H2 
Q{J)=F(i)*F(3)*H2 , 

Q ( JM ) =-F (4 ) *Y( K) *H3-F { 3 ) *U ( K-1 ) *H2+F { 3 ) *U t K) *F ( 2 ) *H2 
20 CONTI 
■ -113 IJ = II-i 
IK=I 1-3 
B( IJ,IK)=0. 

B( IJ,IK+i}=-Hl 



B(ij,iK+ 

RC IJ,] K+.. )=C.3 
Q ( IJ ) = . 

Bdl ,IK)=F(;. )«H'' 

3 U' I » I r>+ ) = F (L ) -'I'll' 

B( II ,lK + „_)=C..j 
oC II , 3 1 )=-h7^:=H,:-h: 

will ) F ( “r ) ^ Y ( j ? j ) sj! }-| J — p ( ) jjcl I / I ^ 

RETUKV ■ ' 

7NC 



as FTC 


1: 



SiJl-R'-'"-. ! JMhN,0,R,X) 

C U.I M (, 1 '‘1 / T ‘ J / U f Y j S I G V f S I G y 

3 ij ’ ' = ■ ^ ^ ^ » S I GV ( 2 ) , D ( 40 , 40 ) , K ( 4 C } , X ( 4 0 ) 

SUf-'XC= ! ” 

SUl-’X 
SUMX J = 

SiJ>!X :ij=>.'. ■■; 

S!JMIJLi= . 

SUMK2J= . 

A X i J A = f. , ■„ 

AX2JX2=0.0 
A X i... J iJ“ • 
fli Y X — • 

H;. =..0/SIGV(:.)'i“|t£ 

Hi=l ,:(/SIG*aI’!''Y2 

oa 1 ■ 1=. ,NN.. 

J = 2 * I “ J. 


K=j+:. 

L = K+ i 

SUMX;. = SUMX.+X( 

SUMX?.= SUMX2 + X(K)’)'*2 
SUMX J, 2 = SUMX 1 2+ X { J ) YX ( K) 
SUhXiU=SUMXlU+X( J)*U(I ) 
SUMXZU = SIJMX2U+X{K)*U( I ) 
SUMUU=SUMUU+U( I )*UtI ) 
SUMXZJ=SU^X2 J+X(K)*X{K) 
AX2JXI=AX2JX1+X{L)>:'X{ J) 
AX2J X2 =AX2 JX2+X ( L ) ( K3 

AX2JlJ=AX2JU + X{L)*U(i ) 
AYX=AYX+XtK)*Y(I ) 
CQNTINUb 
0( 1, 2) =SUYX1«H„ 

D( 1, 2) =SUMX12*H1 
D( 1 , 3) =SUi«XlU*HI 
0(1,4)='/..;. 

D(2, i)=D(l,2) 

D( 2, 2) =SUMX2*Hi 
D(2,3)=SUHX2U=i‘Hl 
0(2,4)=f .r: 
0{2,l)=0(i,3) 

DO, 2) =0(2,3) 

0(3,3) =SUMIL*H1 
D(3,4)=-v.C 
0(4, 1) =0.0 
0(4,2) =0.0 ' 

n=2*NN 



U { j ^ ) 

= { SU?' X2J + 

R(j) =l\ 

X JX:.«Hi 

R{ 1’) =d 

XX JXX«H. 

'■•1 { .■ ) 

X. ... J U H . 

R ( O = { 

A ¥ ( n ) 

RcTUR'j 


2:\!D 



I ) 


ixN) )«h2 




^=ISFTC 


i Nb 

I \'V‘.RT(N,AA,B,X) 

iru: ! 


W 4. W W i ^ 


SLIBR fn) 1 1 _ r- J K 

SUHROUTli. 

DIME'JSin . , 

p ^ ^ ^ ^ ^ ) » e ( 4 0 ) , X ( 4 -. ) , A ( 40 , A I 

DO 199 I=.,N 
I, 'iO )=9 I 1) 

DO 199 J=^,N 
A{ I, J)=ftA(I,J) 

;.99 CONTIM!)',: 

K=: 


KJ=K-1 
DU Z I = , , i\ 

DO 2 J=:,,N 
IFfK.Mc. ,, JGC TO 
PIV0T=a 3S? A (1, j ) ) 


I i =i 
12 = ?. 

GO TO - ■ 

3 DU 5 KK=*,KJ 

IF (I .EQ.IRIKK) )GQ TO 2 
IF{J.EQ.1C(KK) 3G0 TO 2 
5 CONTI MU= 

PIVOT=ABS(A{i, J) ) 

1-1= I 
I2 = J 
GO TU 

2 CONTINUE 

4 DO 6 I=1,N 


DU 6 J=i,N 
IF (K.EO..,, )GC TO 7 
DO S KK=l,KJ 
IF(I .EW.IRCKK) )G0 TO 6 
IF(J.EQ.IC(KK) )G0 TO 6 

8 CONTINUE 

7 IF(PIV0T.GS.A8S(A{I, J)))GO 70 6 
PI VOT=ABS(A{I, J) ) 

Il = I ■ 

I2=J 

6 CONTINUE 
IR(K)=I1 
ICIK)=I2 
IJ-N+l 
DO 9 J = ?.,IJ 

IFU.EQ. I2:)GQ TO .9 ,, 

A(I1,J)=A{I1,J)/A{I1,I2) 

9 CONTINUE 

AI li ,12) =i.C 
DO lu I=:,N 


MATRIX 



D'3 J = ,,,iJ 
I F { I . r Q . 1 R ( l< ) ) jQ 
I F ( J » I: iv . i 1, ( K ) ) GO T Q 

A ( I , J ) = ( I , J ) _ ^ ( j ^ t p 

COOTIO'J. 

O'.) GO I = y ,\ 

I F { I . - 0 . 1 „ ) G ij T f 1 8 9 
4{ I, I? )=l;. 0 
CONTIMUi, 


> 0 

>*^ { 1 1 » J ) 


K=K+ 1 


IF(K.L.::-. J)GC TiJ 

ou :8 K-:.,u 

IL=I R{ K) 

!M=I C{ K ) 

X( IM )= '■,{! L, NO) 
GOrjTIMU;; 

RGTURM 




IBFTC 


SU'’?niJVI 

USING dynamic programming 

^ '■ ^ » S I GM , X EL ECT , S I GE , X C I V I L , S I GC 
T . c I , ! ' c ■' , . ; ^ ‘ ^ ^ ° ^ { 10 ) , s IGE ( i 0 ) , xc I V I L ( 1 0 ) 

D1 ML ' lb i t.'N o IIjL { 10 ) 

CUMMG']/I1TK/C1,C;: 

C ' I M M H \i / 1 i T / 1 S c -L I , I M I T A , i-s N 3 

CiJf^MJN/,.R?A3/!J..,US,ii ,L1,L2,L' 

CJMM':jN/.a'v i3/UtJ 
CUMMON/AK . ^.>/N 
COMMUN/AR'=a5 71MN 

I NTE Gt R IJ ... , U2 , IJG , , X2 , XO , LU_, , UU2 » UU.; , ROUTE t PATH t RATH 

DI MtNSirjiJ PA rH(£,22,2.4) » RATh ( 22, 22 ) , ROUTE ( ,10 ) , Fi ( 22, 22 ) ,F2(22, 

I 77 ) ' ■ 

DiMENSlnr. ISSAT(2),1MITA{2,&),NN3(6) 

IMN=0 


PRINT?.;:.' , (I'SLAT(n,I = i,2) 

25- FORMAT{. X , 1 S -AT* ,21 7 ) 

PR I N 7252 , ( { I :il T.A ( I , J ) , J = l ,6 ) , I =1 , 2 ) 

252 FORM.AT ) 

PRINT2 53, {NK3{ I ),!=;, ,6) 

253 FORMAT {l X,*.\N3Xc, 617) 

N = 1 


IJK=1 

Ul=IMITA{ IJK,1) 

U2 = IMITA{IJK,,2) 

U3 = IMITA(IJK,3) 

L1=I MITA( I JK,4) 

L2=I MITA( I JK,5 ) 

L3=IMITA(IJK,6) 

NNi=Ul-L,:. + l 
NN2~ Li ~ i 
PRINT289 

289 FORMAT(/2 ,:X,«DYNAMIC PROGRAMING*/) 
DO r.i l=i,NN.. 

J = 0 

Xl=i+NN2 

JSEAT = ISEAT(IJK)-X1 


X2=L2 

X3=JSEAT-X.2 ■ ■ -■ 

IF{X3.GT.U3)G0 TO 103 

IF(X3.LT.L3i60 TO 

IF(X2.GT^U2)G^^^ 

JK=X2-L2+i 

j:==u+:r:'^' ' "S ■ 

FI (I ,JK)=FUM(N,X1) + FUN2(N,X2)+FUN3CN,X3) 
IMN=1 ' 

L03 CONTINUt 


x 2 =x:+ . 

GO I;J 
CONTiMij^.: 

I JK= ';:, 

1,3 W = N+* 

t;:=i Mj T.\( i jK,:, } 
iJ2 = IMITA( i JK,.:) 

1J3 = I "'ilTi ( ! JK,3) 

(i JK,4) 

L2 = IMITa{ IJK,5 ) 

L:> = I HI TA ( I .JK,5 ) 

K.JI=i 

IF { UK. cQ. ..)KJ 1 = 2 

301 : 'j N 1 = !,! } ~ L + .■ 

DO j,.2 
NN=H-a 

J=. 

Xi =I+Mri,'. 

JSEAT=ISt;AT{IJK)-x: 

X2=L2 

14 X3=J ScmT-XI 

IF {U3.lt. X3)GO TO .55 

if{l:;.gt.xii)gu to u 

IF(U4.LT.X2)G0 to ^5 
JK=X2-L2 + 1. 

F2 ( I » J K ) = FUN A ( N » X I J +FUM2 ( N, X2 I+FUN3 ( N , X3 ) 
IFttUM: (U-X..).LT.lMITA{KJIf4))G0 TO 16 
LLi = ?'4^j3( D-Xi 
GO TO n 

16 LL1=IMITA{KJI , 4 ) 

17 IF( (NM2{2)-X1).GT.UITA(KJI, ,,))G0 TO 18 

UU1 = MN3{::}-Xi. 

GO TO 19 

18 UU1=IMITA(KJI,1) 

19 IF({NN:Ma)-X2).LT.IMITA{KJI ,5>)G0 TO 20 

LL2=NM3{3)-X2 

GO TO 21 

2r| LL2=IHITA(KJI,3) 

21 IF((NM3(4)-X£) .GT.IMITAtKJI ,2) )G0 TO 22 

UU2=NN3(4)--X2 
GO ITO 23 

22 UU2=1MITA(KJI,2) 

23 IF{(NN3(5)-X3) .LT.IMITa(KJI,6) )G0 TO 24 

LL3=NM3:{5)-X3 , : 

: GO ' TO \\ 24Cr : 

24 LL3=IH1TA(KJI,6) 

240 IF ((\\i3{6)~X3) .GT.IMITA(KJI ,3nG0 TO 

T.-' v>/: UU3 = NN3 ( 6 )-X3 
•; GO TO 26 


25 


25 (Kji,,-.) 

26 I J =0 

no ,?7 k.. = ll:,uu:. 

I I=Kl+.-l;*ITA(!<lJI,.: ) 

JJSF. :T^IS..,.T(KJT )-k 
K2 =L L2 

28 K3=J JSi: /.,l-K2 

I F (K5. GT.LU'J )ijO TT '• 

) T j "h 
IF (K 2, GT.!jb2 )G J ir. 2 ^ 

JKi = K2 + :,-lMlTA(KJI ,8) 

I J = I J+-. 

IFdJ.GT.:, )GQ TO : 

, CrjST=F„(iIrJKi) 

RuUTr ( j =K:.. 

ROUT8d:)=K2 
GO TD v 

30 I F {COS T. LG. F: {II, JK. ) )GU TO 3 
COST=F. (II,JK.) 

ROUTc(:!=K„ 

ROUTh{ 2) =K'-. 

3i K 2 =K 2+ i 

GO TD OS 
29 CONTINJUE- 

27 CONTINUE 

F2(I »JK)=F£ (I, DK)+COST 
PATT(N,I,JK)=RuUTF( ) 

RATH{N,I r JK}=RaUT2(2) 
i55 CONTINUE 
X2=X2+I 
GO TO 14 
15 CONTINUE 
i2 CONTINUE 

DO 40 1=1, HNI 

Xl=I+N!N?. 

JSEAT=ISl;AT{1JK)-X.. 

X2=L2 ■ 

41 X3=JSEAT-X5 : 

IF(X3.GT.U3)&0 TO^ 

IF (X3.LT.t3IG0 TO 44 
IF{X2.GT.U2)G0 TO 40 
JK=X2-L2+. 

FI (I ,JK)=F2(I» JK) 

409 CONTINUg 
'X2=X2+3 

, GO " TO : 41, ..'v O, 

40 CONTINUE 

IF(N.EQ.5)GC TO -Oil 
IFd JK.FQ.DGO TO 304 



GiJ T1 5 
3 ' 4 I J K— i 
305 CONTIMUc 

Gtj Tj ^.3 
3 i.. CO.MI'ilL 
I J = 0 

DO 5C i = 

JSl:Al=IS..-iT{IJK)-X.' 

X2=L2 

51 X3 = J S - 1 ~ A ‘i 

I F (X3. GT.Li )(;o TO 53 
I F ( X :5 • L T • L it ) =j Li 1 U 5 
I F { X • kj r • U ) GO To 50 
JK=X2-L1 + ., 

IJ=IJ+„ 

IF(IJ.GT..)Gn TO 52 
CaST=FI ( 1 , ,JK ) 

'TauTc{.:.)=x:. 

RUUT3{i;)=X> 

GO TO 53 

52 IF(COST.Lc.Ft ( i, JK) )G0 TO 53 
cost=f:; (I, jK) 

ROUT .:;(!) =X1 
H0IJTe(i:)=X2 
53 X2=X2+; 

GO TO ■ SI ' 

50 CONTINUE 

PRIMTii.. ,C0ST,R0UT£(1 ),R0UTE{2) 

iOO FORMAT (//30X,^CPT. COST = »!=, F 15. 8, 5X ,*ROUTE ( 1)=*, I3» 5X, #ROUTE ( 2 ) =* » 
ii3//) 

Ml=ROUTE(i) 

M2=RUUTE(;:) 

I=M1-Ll+UN1 
MN=M2-U1+NNI 
M3=PATHt5,I ,M:i) 

M4=RATH(5,I,MM) 

I=M3-L2+NNi 
MN=M4-U2+WN1 
M5=P ATH(4, 1 
M6=RATH(-.-il fMN) 

I=M5-Ll + NN1 
,MN=M6-Ul + NN1 , 

M7=PATH(3,I ,MN)t 
M8=RATH( j ,I ,MN) 

: 0 .i=h7-lo+n'n: t 
MN=MS-U2+I\N1 
^ • H9=PATH(2,I tMM,) 

vt/O M10=RATH{2,I,MN) VO-:: v : : 't. ; 

PR I M T1 C* j M I , M£ I M3 , M4 » M5 1 M6 , F7 j M8 , M9, MiO 


iO i 


F11Rv,^\T(5 

Fr^lD 


R -SUITS*// :0(5X,I 5 )///5CX, ^PROBLEM 


IS 


OVER*) 



«IBFTC 

FUi'iC 1 I ‘ - iv; F L', i!'.. ! *'] » U!'» ) 

C COMPOSIT:-^ TRAPEZCIL RULE IS 

{, iJ ..VuLiJi; rii;;. FUfCTION VALUE 

C FUNi 
C 

COMMn"J/I ilK/'- ,C2 
crjyMu^j / A R., A 2 /u . , UE , u-: , l i , l2 ,l 
I i'jTi G” K j L Z j U3 » U:\'» UNI t UUI 

Fx=F a. ,UN) 

IF (L U.U:) )G'J TO 30 
L L .1. = L j + ^ 

DO i i~LLijl I'j 
Ak Fl=F:i+:., *FK, (I, UN) 

3C CONT I MUl: 

Fi = ".5«F.. 

F2=FN1 (U.,UN) 

IF (Uj . LQ.U?0 GO TO a- 0 
UUl=Ui-i 
DO E' I=UN,LU1 
20 F2=F2+0. «FNZ.{I,UN) 

4C continue 

F2=-F2/2.C 
FUNl=:;: *Fi, + CA*F2 
RETURN 
END 


USED 


TO 




❖IBFTC 


FUMCTI FL' i.; ( 1,UN) 

COMHOM/I ITK/C. 

C 3 M n \i / A ■< .- i 2 / l) f U » U , L i , L 2 , L 

I ''JTl-G-*'’ G »L -»U’..,Uil,UNl,UUi »UW3,UU3» UU2| UN2 
Fi =F A? ( L . , Uh ) 

IF ( L f'* • U ) Gu To 

DO 1 i=LL:,lJO 
^0 Fi.=F;' + 2 .1 ^'FM( I ,UN) 

3i.. CONTIMLih 
F 1 =^F 1/ " . 

F2=FH2 (U. jUN! 

IF(U?.LQ.L'nGO TO 40 
UU2 = U2-,. 

00 2' I=U\,l!L 
2:' F2=F2+2. *fn:: ( I ,UN ) 

4u CONTINUE 
F2=-F2/2. : 

FUN2=Ci«Fl + C:'t*F2 

RETURN 

FNC 



FUMCTI L:.H Ffi.:' ( M, UN) 
COMMaN/1 ITK/C" jf-;’ 

C0FM0N/A>iA2/U.,.,'u-,UX.,L 

INTEGER U: ,U.;, LJ::>,Uri,UNl 
F1=FN? ) 

I F (13. cO.UK )G!J TO 
LL3=L2+I 
no 3 0 I =113, ON 
Fi=Fl+,.', ^'F^';,■{ I ,U!\ ) 

3 CONTINUE 


E»L2 ,LE 
,UUI,UN3,UU3 


F*=') .5*F 


F2 = F InJB ( U.'i j UN ) 

I F (U 3. EU. UN )G0 TO A 
UU3=U3-1 


DU 2 0 I=UN,L'U3 
£.0 F2=F2 + ;':.0*FN3{ i ,UN) 

4. CONTINUE 
F2=-i . 54F2 
FUN3 =C.,.«F^ + C2^F'. 
RETURN 
END 



>i*IBFTC 


FUf<:Tib;4 F^.JXTUM) 

C0MM3'UAREa5 ■/ IHN 

COf^^'iO^/ARiAl/XM, SlGiM,XE,SIGE, AC, SIGC 
C 0 M M 0 M / A R ,4 2 / U , U ^ , U 3 , L 1 , L 2 , L i 
COMMON /AR3AS 12/ ^ 

DIMCFslSlUO XOECH(IQ) jSIGMdC) ,XELECT { 10 ) , S IGE { 1 0 ) , XC I VI L ( 10 ) 
DIMENSlC'j SIGC (1C) 

0 1 Ml 0 S 1 0 'i X M ( ) , X 3 { i ), XCl I ) 

I N TE Gr R U N , X , XH rl CH , X S L ECT ,X C I V I L 
DO aoo 1=1,5 
XMECH( I ) =X'-* ( ! ) 

COF'T I N IJl 
iNj = e-M 
Ci=U'N-X 
C 2 =X-xr,:CH{N) 

F1N1 = C1 *3XP ( - : . 5*C2^*3/S IGM ( N) ) 

RETURN 

END 













^IBFTC 


FUNCTI FM! ( X, UN ) 

CaMMOM/AR^45 /IMU 

CO.yMON / /.R S A 1 /X M , S I GH , XE , S I G E , XC , S I GC 
!j I Me :^iS I u:x X ( : ■) ) , Xc { ^0 ) , XC ( ; :) ) 
COMMiliN/AX .A2/U...,U:,.',U3,L1,LZ ,L-> 
COMMOM/ARXAlZ/i'l 

DlME'slSIU ! i . ) , SIGM{ IC) ,XcLECT{ 

UiMEUSIU*) S!Gt(iO) 

I N T & ■ R UN, X , X M s CR , XEL ECT ,X C I V I L 

DO 3:' I = :i,5 
XELcCTd )=XE{I } 

CONTINUE 
N=e-M 
Cl=UM-X 
C2=X-X. Led (M) 

FHL=C;i*EXP{-u. 5*C2=«=«L/SlGt ( M I 
RETURN! 

END 


iC),siGEtia),xciviL{ia) 



*IBFTC 


FUfJCTIfiN F^ 3 ^,(tUN) 

COMMOi'J/ARcASo/ 

COMMON^ /.\K ; SI GM , XF , S IGfc 'C cirr 
DIKE.'vSU.-, Mil. 

DIME'^SIOa SIGC(l“i 

li'iTcGirR B'N » A » XH£ CH f X*"- L ’-fT y'~rv/Ti 

DO 300 I = i,5 

XCIVIL (I )=XC(I ) 

CONTIMUF 
N = 6- M 
Ci=UN-X 

C 2 =X-XCI VILIM) 

FN3=Cl>i'£XP(- .5=}:r2:?=* VSIGCfM ) 

RETURN ^ 

END 


) f XCIVILI 1 C 



❖IBFTC 


SUBRULlTi !j£ AV -R AG { X , XX , Nl i fFF,SIG) 

COr'^MON/TAJ/U,Y,SIGV,SIGW 

Ca.YMDN/AREAl:;/IT 

DI KEf'lG 10.4 X { ■'i.i'.) ) f X >;{ ' 0 ) , SI G (iO ) , FFF( 12 j4 ) 
DIMENSION I.U(5),U(40),Y{40),SIGV(?) 
RFADl" , (UU{ I ) , I=„ ,5 ) 

F0RMAT(5F1 : .J) 

N = 14 


IF (IT. .)N = :.6 

NNi = '-.* (M, +:, ) 

NN=2’:'N! 

MH=1 i, 

IF(IT.i.u.,0)YN=„3 

MM=2*MN 

mnN=MN-2 

I F ( N Hi . 0 Q , M ^' ) G a TO 20 
DO IS I=N1,MNN 

IJ=2X=I+£ 


N2=I+l 

IK=IJ-. 

IL=IJ + * 

II=IJ + 2 
X(IL)=X(IJ) 

X ( in = FFF(Nl,i )*X( IK)+FFF(Ni,2)XcX{ IJ )+FFF(Ni,3 )X=U(N2) 
iO CONTINUE 
2i> CONTINUE 
MK=MN-i 

DO 40 I=HK,N 

IJ=2*I+2 

J=I-MN + 2 

IK=IJ-1 

IL = IJ-»-l 

n=ij+2 

XdL ) = X(IJ) 

X( 11 ) = FFF(Nl»i )*X( IK)+FFF(N1»2)*X(IJ)+FFF(N1,3)=«'UU( J) 
40 CONTINUE 

DO 45 1=1,5 

J = 2=«'(I+MNI 
45 XX (I )=X{ J) 

PRINTl 01, (XXd ), 1 = 1,5) 

101 FORMAT (luX, LVALUE OF XX*,5F20.5) 

CALL C0EFF(X,FFF,N1,SIG) 

ReTURN 
■ END 


=S=IBFTC 


COt:FF(X,FFF,Ni,SIG) 

Jt I t , 4 „ 


_SUB^aUTI'.. 

2) ,Q(4,2),5lG(.. ),AB(2,2) 


REAL K 
DO 12 I=i,2 

DO lit J = l,.l 
P{ I, J)=- . 

4{ 1, i) =0.0 
ACltc) = ’ 

A(2,3 )=FFF(M,*) 
A(2,2)=FFF(N:. 

DO 20 1=1,2 

DO 3^ J = l,2 


3) 8(i,J)=L.:: 

2v P(I,I)=1.. 

C(i)=0.0 
C(2) =FFF (N >. ,4) 

D = l.. 

K(l,l)=i,C 
K{2,2)=*2.C 
K{1,2) = , 

K12,I)=' . 

L = 5.0 
IJ=0 

9': IJ = IJ+i 

Hi=C (2 )**.:’«F(2 ,'’)+3. 

H2=1.0/H... 

F{].)=C{2)^P(^,1) 

F(1)=F{1)*H2 

F{2)=FFF(Ni ,i) *C(2 )*Pn ,2) + FFF(Nl,2)*C(2)*P(2,2) 
F(2)=F{2)*H2 
DO 40 i=i,2 

ABd ,1 )=A(I,a) 

40 A{ I, 2)=A( I ,2)-F( n^C{£) 

DO 50 1=1,2 

DO 50' J = X,C 
Q( I, J)=; . 

00 6& KK=1,2 

60 Q( I, J) =Q(I , J)+P(I,KK)*A8(KK,J) 

50 CONTINUE 

DO 7u 1=1,2 
:00 70 J=l,2 
R( I, j) =0.0 
00- 79 KK=1,2 ■ 

7C R(I, J)=R(I,J)+AB(I,KK)*Q{KK,JJ 

' ^ ' : a:3 

■ 1 00-0 BOO d=l;:,2- , 

80 P(I, J>=R(I, J) + K(I,J)+5.0*F{ 

PRIMT aCh^i ( (P( I , J) ,J=lf 2),I=l,21 


F ( 2 3 , 


FORMAT {IX,*F{ I » J ) 2F1 5 ) 

I F{ IJ. LE.i . )GQ TO 9C 
IK = IJ-... 

SIG( IK )=P( 2,2} 

IF{I J« LT»„3 )G(J TO 9| 

RETURN 

ENO 



INPUT MTA 


1. ffiCEANIGAL 

ENGI'TEEPIrlG DTOOTPTtijtt 




\ I 195 

98 

105 74 

134 

56 

120 

148 

71 

119 10 




\ • 50 

25 

50 25 100 50 

125 

60 

100 

50 

100 

50 




Variances: 5,0 

3.0 12.0 




2 » EIEGTEICAL EHGIiTSBRTNO DISCIPLINE ' 




Yk ; 142 

77 

150 88 

157 

80 

136 

70 

120 

65 




\ • ^20 

65 

100 40 60 

30 

65 

35 

60 

30 





Variances: 5, 

0 

3.0 12.0 




3, CIVIL ENGINEEHING DISCIPLINE ' 




\ 

77 

78 68 171 

37 

60 

30 

78 

38 





: 75 

40 

60 30 50 

25 

60 

30 


55 25 

®nax ^®chanical, Electrical and Civil Ei^ineering 
disciplines : 90, llC, 80 respectively. 

®i 3 in Meclianical, Electrical and Civi 1 Engineering 
disciplines : 70, 90, 60 respectively. ' 
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